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ABSTRACT 


A  STUDY  OF  SPECIFIC  METHODS  OF 
MEASURING  LUNAR  GRAVITY 


The  measurement  01  the  gravity  field  of  the  moon  from  an  orbiting 
satellite  is  examined  from  the  points  of  view  of  requirements  for  geodetic 
and  geophysical  application,  and  physical  phenomena  and  instrumentation  to 
meet  these  requirements. 

,  It  is  shown  that  it  is  feasible  to  determine  the  gravity  field  by 
measuring  the  gravity  gradient.  Mathematical  relations  useful  for  data 
processing  are  stated.  Examination  of  several  principles  for  sensing 
gravity  gradient  leads  to  the  recommendation  of  the  vibrating  string  gradio- 
meter  as  the  best  implementation  with  size,  accuracy,  range, and  versatility 
being  in  its  favor.  In  addition,  the  engineering  status  allows  prototype 
design  to  begin  immediately. 
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SECTION  1 
INTRODUCTION 


The  purpose  of  this  program  is  to  study  the  feasibility  of  measuring 
the  gravity  field  of  the  moon  from  an  orbiting  satellite.  The  utility  of  such 
data  lies  in  their  application  to  geodesy  (determining  the  shape  of  the  moon) 
and  to  geophysics  for  a  better  understanding  of  the  structure,  density  and 
internal  processes  of  the  moon.  In  immediate  application  it  offers  additional 
data  for  trajectory  calculations  for  soft  landing  and  manned  exploration  of 
the  moon.  The  principles  are  also  applicable  to  better  definition  of  the  earth 
gravity  field  for  precision  navigation  of  missiles,  and  for  scientific  planetary 
probes. 


Direct  methods  for  obtaining  these  data  require  landing  a  gravity 
meter  on  the  surface  of  the  moon  and  a  geodetic  reference  to  relate  points  of 
gravity  observation.  Without  mobility  the  utility  of  even  multiple  meter 
landings  is  severely  limited,  but  the  difficulties  of  putting  a  vehicle  on  the 
moon's  surface  are  too  well  known  for  further  discussion  of  this  as  a  technique. 
At  some  future  date  this  will  be  both  desirable  and  necessary. 

The  indirect  methods  involve  analysis  of  the  trajectory  of  an  orbit¬ 
ing  vehicle  or  the  measurement  of  gravity  gradient  on  board  the  vehicle. 
Investigation  of  this  latter  forms  the  basis  for  this  report.  Theoretically  the 
measurement  is  sufficient  to  allow  description  of  the  gravity  field,  and  the 
environment  of  the  measurement  ic  a  low-noise  one  suitable  for  the  measure¬ 
ment. 


Instrumental  design  analysis  further  indicates  it  is  feasible  to  make 
the  necessary  measurement.  A  vibrating- string  gravity  gradiometer  is 
presented,  which  conservative  estimates  indicate  will  allow  anomaly  deter¬ 
mination  to  the  10  milligal  level;  if  the  instrument  preserves  the  dynamic 
range  of  its  predecessor,  the  Arma  vibrating  string  accelerometer,  it  will 
allow  anomaly  mapping  at  the  0.  1  to  1  milligal  level.  Because  this  mapping 
is  global  and  fast  it  appears  to  offer  an  ideal  solution  to  lunar  gravity  deter¬ 
mination. 
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SECTION  2 


CONCLUSIONS  AND  RECOMMENDATIONS 
2.  1  Conclusions 

Using  gravity  gradient  measurements  it  is  feasible  to  determine 
the  gravity  field  of  the  moon  from  an  orbiting  satellite.  The  level  of  per¬ 
formance  with  which  this  can  be  accomplished  provides  the  following  informa¬ 
tion: 

a)  sufficient  data  to  permit  definition  of  the  shape  of  the  planet 

b)  at  least  two  orders  of  magnitude  improvement  in  the 
accuracy  of  current  knowledge  of  the  strength  of  the  field 

c)  mapping  of  anomalous  structure  on  the  selenoid  at  the  one 
to  ten  milligal  level 

It  has  been  assumed  that  the  probable  orbit  will  be  about  100 
miles  above  the  lunar  surface.  In  recent  discussions  there  have  been  in¬ 
dications  that  a  forty  mile  orbit  may  be  employed,  in  which  case  the  anomaly 
discrimination  level  will  be  further  improved  by  more  than  tenfold.  Apart 
from  generally  improving  with  reduced  orbital  altitude,  the  system  does  not 
constrain  the  orbit  in  any  way.  Indeed  we  note  that  the  natural  eccentricity 
and  precession  of  the  orbit  provides  the  global  covering  sample  most  useful 
for  evaluation. 

The  data  processing  problem,  the  computation  of  gravity  field 
from  the  gradient,  is  amenable  to  ordinary  computer  technology.  Some 
computational  techniques  are  suggested  to  illustrate  feasibility  in  the  report. 

As  primary  instrumentation  the  vibrating  string  gradiometer  is 
a  more  than  satisfactory  instrument  for  this  purpose.  A  preliminary  design 
analysis  and  performance  summary  is  presented;  the  performance  is  predi¬ 
cated  upon  the  Arma  vibrating  string  accelerometer  and  the  design  is  an 
engineering  extrapolation  of  the  knowledge  and  experience  gained  from  the 
low  range  accelerometer  now  under  test  for  RTD.  A  minimum  sensitivity 
(dynamic  range)  of  one  part  in  1 0 5  or  better  is  indicated,  compatible  with 
the  anomaly  mapping  goals. 

As  designed,  the  instrument  also  incorporates  a  low  range 
accelerometer  capability  which  may  be  employed  in  calibration  and  for 
ion-thrust  engine,  drag  and  solar  radiation  monitoring. 
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Note  that  the  instrument  is  readily  calibrated  in  orbit  with  high 
precision.  Biases  in  the  instrument  are  automatically  isolated  in  the  data 
reduction  process,  including  time  dependent  drifts  of  the  bias.  The  scale 
factor  is  readily  determined  by  several  methods,  the  simplest  being  the  use 
of  the  gravity  field  of  a  nearby  positionally  controlled  mass. 

2. 2  Recommendations 

In  accord  with  the  above  conclusions  the  following  is  recommended: 

1  -  Detailed  design  and  fabrication  of  a  vibrating  string  gradio- 
meter  should  be  carried  out.  The  gradiometer  should  be  fabricated  for 
satellite  application.  The  test  program  initially  should  include  zero-g  air¬ 
craft  flight  and  recoverable  flight  testing  in  sounding  rockets.  Particular 
attention  should  be  given  to  the  possibility  of  designing  a  proof  mass  suspen¬ 
sion  system  which  would  allow  fully  sensitive  operation  of  the  gradiometer  in 
a  one-g  cross  axis  field.  If  this  is  achieved  the  gradiometer  can  be  employed 
on  a  horizontal  platform  to  measure  gravity  anomalies  of  the  earth  from  an 
aircraft  in  level  flight. 

2  -  A  complete  computational  system  should  be  worked  out.  In 
particular  numerical  methods  for  calculating  the  anomaly  field  more  directly 
from  gradient  measurements  than  those  indicated  in  this  report  should  be 
developed.  These  will  be  of  use  in  satellite,  airborne  and  ground  based 
gravity  surveys. 

3  -  Further  detailed  study  of  the  operational  usage  of  the  gradio¬ 
meter  and  the  system  should  be  made.  This  can  be  done  partly  by  simula¬ 
tion  to  examine  alternative  methods. 

4  -  The  gradiometer  provides  information  on  the  gravity  field  in 
different  forms.  For  example,  the  curvature  of  the  field  is  directly  available. 
It  appears  desirable  to  explore  direct  methods  of  interpreting  these  findings 
as  an  alternative  to  the  conventional  anomaly  maps.  This  entails  deeper 
investigation  into  the  geophysical  interpretation  of  the  data. 

5  -  A  study  should  be  made  of  the  gravity  gradient  technique  as 
a  possible  experiment,  and  as  a  working  tool  on  a  manned  orbital  research 
laboratory.  The  presence  of  the  human  monitor  allows  a  number  of  opera¬ 
tions  to  be  carried  out  semi -automatically ,  and  data  can  be  evaluated  on  the 
spot .  If  the  satellite  is  properly  equipped  the  gradiometer  could  even  be 
assembled  in  space  and  tension  and  proof  mass  variations  could  be  tried. 
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The  use  of  the  gr&diometer  to  indicate  the  local  vertical  can 
facilitate  active  attitude  stabilization  of  the  satellite  laboratory.  At  the 
sensitivity  level  of  the  proposed  instrument,  and  with  the  aid  of  the  accel¬ 
erometer  function,  short  range  non- radiating  docking  control  functions  can 
be  carried  out.  Drag  measurements  can  also  be  made,  and,  of  course,  the 
fundamental  geophysical  data  that  gravity  offers  can  be  obtained  and  inter¬ 
preted. 
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SECTION  3 

TECHNICAL  DISCUSSION 


Introduction 


This  section  contains  the  important  analytical  and  design  findings 
of  the  study.  They  are  presented  in  the  following  sequence: 

a)  The  use  of  gravity  data  is  first  discussed.  The  magnitudes 
of  gravity  anomalies  on  the  earth  and  typical  geophysical 
phenomena  are  cited.  The  current  status  of  knowledge  of  the 
moon  and  other  planetary  bodies  is  then  summarized.  Sen¬ 
sitivity  requirements  for  gravity  data  to  be  useful  are 
determined. 

b)  Equations  of  the  gravity  field  are  derived;  relations  to 
measurements  by  gradiometer  are  derived  and  methods  are 
presented  for  deriving  anomaly  maps  from  the  observations. 

c)  Instrumentation  techniques  for  measuring  the  gravity 
gradient  are  reviewed,  and  a  detailed  analysis  of  a  vibrating 
string  gradiometer  with  preliminary  design  data  is  shown. 
Operational  configurations  for  calibration  and  orientation 
are  discussed.  These  are  related  to  the  requirements. 

3.  1  The  Dependence  of  Geodesy,  and  Geophysics  on  Gravity  Data 

The  gravity  field  of  a  planet  is  a  unique  function  of  the  distribution  of 
mass  in  the  planet;  that  is  it  depends  upon  the  shape,  density  and  position 
of  matter  in  the  body.  Through  the  process  of  hypothesizing  density  distri¬ 
butions  and  mass  transfer  mechanisms  and  then  testing  the  predicted  gravity 
field  against  the  observed,  the  geophysicist  is  able  to  determine  the  most 
probable  internal  processes  of  the  planet. 

The  uniqueness  of  the  gravity  field  has  been  fundamental  in  the  devel¬ 
opment  of  navigation  techniques.  The  position  of  the  local  vertical  projected 
on  to  the  celestial  sphere,  and  the  time  of  day,  suffice  to  determine  location 
on  the  earth.  To  effect  this,  of  course,  the  navigator  must  know  the  gravity 
field  everywhere;  the  standard  representation  of  this  is  a  map  -  a  correla¬ 
tion  of  geographic  detail  with  coordinates  on  the  celestial  sphere.  However, 
in  order  to  determine  lineal  distance  between  points  on  the  map,  one  must 
put  the  map  not  on  the  celestial  sphere  but  on  the  real  solid  figure  of  the 
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earth;  the  geodesy  of  the  planet  must  be  established.  In  determining  the 
shape  of  the  planet  gravity  again  plays  an  important  role.  Strong  correspon¬ 
dence  between  the  geoid  and  the  gravity  equi potentials  leads  to  a  relatively 
simple  description.  For  the  earth,  mean  sea  level  is  the  equipotential  and 
corresponds  reasonably  well  to  the  ellipsoid  now  generally  employed  for  the 
shape  of  the  earth. 

Nevertheless,  current  maps  are  imprecise.  As  Heiskanen  points 
out  in  the  Handbook  of  Geophysics  the  incompleteness  of  existing  triangu¬ 
lation  surveys  leads  to  many  geodetic  systems,  depending  upon  where  the 
geodesist  begins  his  survey.  Since  the  continental  masses  are  not  now  well 
tied  together,  the  survey  can  not  be  closed  globally.  The  data  are  good 
enough  for  ships  but  not  as  good  as  desired  for  highly  precise  navigation  for, 
say,  missiles. 

Accordingly,  wide-ranging  gravity  data,  extended  over  the  oceans, 
and  into  otherwise  inaccessible  regions  has  significant  military  and  com- 
merical  value  to  the  United  States.  Similar  considerations  apply  to  the  moon. 

In  addition  to  its  value  for  navigation,  gravity  data  is  of  fundamental 
interest  to  the  geophysicist.  Basic  theories  of  planetary  origin  and  struc¬ 
ture  can  be  tested  by  observing  the  gravity  fields.  Keys  to  local  sub¬ 
terranean  geology  of  commercial  importance  are  provided  by  gravity  surveys. 

In  this  program  we  are  interested  in  both  lunar  geodesy  and  lunar 
geophysics.  The  following  section  discusses  two  of  the  geophysical  aspects 
of  the  earth  in  more  detail,  as  possible  cues  to  the  geophysics  of  the  moon. 
These  are  the  concept  of  isostasy  and  the  nature  of  anomalies. 

3.  1.  1  The  Concept  of  Isostasy 

During  the  course  of  making  precise  measurements  it  was 
noticed  that  the  attraction  of  large  masses  such  as  mountains  was  less  than 
predicted  by  the  law  of  gravitation.  It  was  concluded  that  the  density  of 
the  mountains  is  less  than  that  of  the  surrounding  rock  (this  is  still  the 
present  day  conviction). 

This,  in  turn,  led  to  the  concept  of  isostasy  in  which  the  earth's 
crust  is  considered  to  "float"  on  a  liquid  interior.  The  mountains,  lowlands, 
and  oceans  are  thought  to  exert  the  same  pressure  at  depths  not  far  (about 
114  kilometers  *)  below  sea  level.  These  masses  are  said  to  be  in  "isostatic 
equilibrium"  and  the  depth  at  which  the  pressure  is  considered  to  be  uniform 

is  called  the  "depth  of  isostatic  compensation.  " 

• 

*  There  are  other  estimates  of  this  level.  The  figure  cited  above  is  due  to 
Hayford.  Heiskanen  has  proposed  a  level  ranging  from  30  to  40  km,  varying 
over  the  earth. 
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There  are  several  isostatic  systems,  differing  in  the  way  in 
which  it  is  assumed  that  mass  compensation  is  achieved. 

In  one  system  it  is  assumed  that  mountains  have  low  density 
earth  crust  directly  beneath  them  and  that  oceans  have  high  density  earth 
crust  directly  beneath  them  in  such  a  way  that  the  higher  the  mountain  the 
lower  the  density  of  the  earth's  crust  and  the  deeper  the  ocean  the  higher 
the  density  of  the  earth's  crust.  If  the  earth's  crust  is  assumed  to  be  sub¬ 
divided  into  columns  floating  upon  the  molten  inner  mass,  then,  according 
to  this  assumption,  all  columns  at  the  same  latitude  have  the  same  total 
mass.  Consequently,  the  bases  of  all  columns  are  at  the  same  radial 
distance  from  the  center  of  the  earth,  resulting  in  a  uniform  depth  of  iso¬ 
static  compensation. 

In  another  system  it  is  assumed  that  the  earth's  crust  is  sub¬ 
divided  into  columns  as  before;  the  columns  now  have  unequal  masses  and 
therefore  float  upon  the  molten  inner  mass  at  different  depths  of  displace¬ 
ment.  Thus  in  this  system,  the  depth  of  isostatic  compensation  is  greatest 
under  mountains  and  least  under  oceans. 

In  still  another  system  it  is  assumed  that  the  strength  of  the 
earth's  crust  is  such  as  to  prevent  a  floating  system  from  being  in  complete 
equilibrium.  This  will  result  in  the  large  mass  of  a  mountain  being  com¬ 
pensated  for  in  a  region  greater  than  its  direct  vertical  projection.  The 
crust  below  sea  level  is  assumed  to  behave  somewhat  like  an  elastic  plate 
floating  upon  a  liquid  and  loaded  locally  by  mountains,  valleys,  and  oceans 
a  negative  load. 

3.1.2  The  Concept  of  Anomaly 

If  the  earth  were  considered  to  be  a  perfect  oblate  spheroid 
without  nonhomogeneities,  all  equigraviational  surfaces  would  be  concentric 
with  the  sea  level  datum.  It  would  then  be  possible  to  calculate  gravity 
very  simply  everywhere.  However,  as  has  been  noted,  the  measurements 
of  "g"  have  differed  by  varying  amounts  from  the  values  calculated  for  an 
idealized  earth  as  described  above.  These  differences  are  attributed  to  the 
nonhomogeneous  density  distribution  of  the  earth's  crust  and  are  called 
anomalies.  The  concept  of  isostasy  was  created  in  order  to  explain  these 
anomalies,  and  the  various  isostatic  systems  described  have  been  used  one 
at  a  time  or  in  combination  in  order  to  minimize  them. 

When  elevation,  topographic  and  isostatic  corrections  have  been 
included  in  a  particular  theoretical  determination  of  "g"  on  the  earth's  sur¬ 
face,  the  residual  anomalies,  determined  by  comparison  with  actual  measure¬ 
ments,  are  relatively  small.  In  the  United  States  they  are  less  than  ±50 
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milligals  (1  gal  =  1  cm/sec^)  with  few  exceptions,  the  average  determination 
being  much  lower.  This  corresponds,  roughly,  to  about  ±50  parts  per 
million  of  Mg"  demonstrating  that  the  earth  exhibits  a  high  degree  of  isos¬ 
tatic  adjustment.  In  the  few  exceptions  mentioned,  the  magnitudes  of  the 
anomalies  may  be  much  higher,  running  as  high  as  hundreds  of  milligals. 

It  has  been  demonstrated  that  some  of  these  occur  in  regions  of  abnormal 
"tectonic  activity.  Presumably,  this  would  indicate  a  temporary  deviation 
fr6m  isostatic  equilibrium. 

Typical  residual  gravitational  anomalies  on  earth  are: 

Complete  Tidal  Cycle 
Small  Faults  or  Ore  Bodies 
Salt  Domes  or  Buried  Ridges  - 
Small  Craters 
Mountain  Ranges,  Ocean 
Deeps,  or  Rift  Valleys 

Gravitational  anomalies  as  determined  from  measurements 
aboard  an  orbiting  vehicle  do  not  uniquely  fix  the  shape  of  the  terrain  below, 
since  a  mountain,  for  instance,  and  a  subsurface  high  density  mass  can  both 
produce  positive  anomalies  of  the  same  magnitude.  On  earth,  topological, 
acoustical,  magnetic,  tectonic,  and  radiological  data  are  used  to  minimize 
the  residuals.  These  have  led  to  conclusions  regarding  crustal  and  sub- 
crustal  phenomena  such  as  stress  releases  in  the  crust,  tilted  fault  planes 
(steps  in  the  crust  caused  by  shifting  earth  masses),  and  broken  formations 
(separations  in  the  crust  caused  by  converging  fault  planes).  If  it  is  true, 
as  has  been  conjectured,  that  the  moon  is  a  cold  body,  that  is  with  a  solid 
interior,  any  evidence  of  crustal  and  subcrustal  phenomena  as  we  know  them 
on  earth,  will  be  useful  in  drawing  conclusions  regarding  the  origin  of  the 
moon.  Under  these  circumstances  there  would  be  no  need  for  concern  re¬ 
garding  the  constancy  of  geodetic  determinations  over  a  period  of  years,  since 
the  concept  of  isostasy  would  be  inapplicable. 

3.1.3  Gravity  Measurements  on  the  Moon 

This  study  is  concerned  with  the  feasibility  of  determining  the 
gravity  field  of  the  moon  from  an  orbiting  satellite.  Before  going  onto  the 
satellite  a  few  words  about  other  methods  is  in  order. 

The  simplest  concept  is  to  land  a  gravity  meter  on  the  surface 
of  the  moon.  This  is  technically  feasible.  However,  the  data  would  have 
limited  utility  since  it  provides  an  isolated  value  of  the  field.  If  properly 


0.  2  to  0.3  milligal 
0.  1  to  0.5  milligal 
1  to  2  milligals 
Up  to  15  milligals 

Several  Hundred  milligals 
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placed  and  if  sufficiently  sensitive  it  might  yield  some  information  on 
tidal  effects,  and,  from  this,  inferences  can  be  made  about  rigidity  of 
the  crust.  An  improvement  would  be  afforded  if  several  meters  were 
landed  on  the  moon.  But  to  use  the  data  for  geodesy  a  triangulation 
net  would  have  to  be  set  up. 

A  more  meaningful  experiment  employs  a  probe  on  a  lunar 
surface  vehicle.  This  would  permit  data -taking  over  a  wide  area  and  through 
an  on-board  navigator,  a  suitable  triangulation  or  geodetic  net  can  be  es¬ 
tablished.  Currently,  considering  the  rugged  nature  of  the  moon's  topo¬ 
graphy  and  our  uncertainty  with  respect  to  its  surface  composition  (dust?  ), 
such  a  program  appears  out  of  the  question.  Someday  it  must  be  attempted  . 
A  surface  program  will  be  necessary  to  provide  the  geological  detail  re¬ 
quired  for  fuller  understanding  of  the  moon. 

However,  it  appears  quite  feasible  to  place  a  satellite  in  orbit 
around  the  moon  and  telemeter  its  data  back  to  earth.  Accordingly,  the 
central  question  is  the  feasibility  of  instrumenting  it.  for  the  collection  of 
gravity  data.  Two  general  approaches  may  be  followed.  The  first  is  ob¬ 
servation  of  the  trajectory  and  the  determination  of  the  planetary  constants 
by  calculation  of  the  accelerations  that  appear  to  act  on  the  satellite. 

Although  feasible  in  principle,  the  data  processing  involves  the 
differentiation  of  data  that  are  intrinsically  noisy.  Good  observations  of  the 
principal  term  of  the  gravity  field,  that  is  the  major  attraction  term,  can 
be  obtained  in  this  way.  However  this  is  already  fairly  well  known.  It 
appears  highly  dubious  that  anomalies  can  be  effectively  detected.  The 
other  approach  is  that  taken  in  the  remainder  of  this  report;  the  use  of  an 
on-board  gravity  gradiometer  to  provide  the  gravity  data.  The  motion  of 
the  satellite  provides  an  essentially  gravity-free  environment,  and  one  free 
of  environmental  noise  to  a  very  high  degree.  Additionally,  because  of  its 
nature  a  gradiometer  is  not  affected  by  external  forces  on  the  satellite 
(such  as  radiation  pressure) . 

Now  it  is  evident  that  a  very  sensitive  accelerometer  mounted 
in  the  satellite  would  indicate  the  difference  in  gravity  due  to  its  separation 
from  the  center  of  mass  of  the  satellite.  However,  this  point  is  not  easily 
found  and  is  subject  to  change.  Accordingly,  the  gradiometer  technology 
considered  inthis  report  leans  toward  measurements  that  do  not  depend  upon 
knowledge  of  the  c.g.  As  will  be  seen  ,  this  is  quite  feasible. 
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3.2  Gravity  Gradient  Theory 


We  define  gravity  gradient  as  follows.  We  first  find  the  vector 
difference  of  the  gravity  forces  on  two  nearby  points  in  the  gravity  field. 
Then  we  take  the  component  of  this  vector  along  the  line  joining  the  two 
points.  Finally  we  consider  the  limiting  process  which  results  when  the 
separation  between  the  points  passes  to  zero.  Mathematically  we  may 
formalize  the  process  by  writing 


d  r 


3.2-1 


The  result  is  a  scalar  function  of  five  variables,  the  three  variables 
of  position  and  two  variables  defining  the  direction  along  which  the  gradient 
is  taken. 


In  this  section  we  shall  start  with  the  gravity  potential  function  V  and 
set  down  the  important  relations  between  V,  the  vector  gravity  (g),  and  the 
gravity  gradient  (gg). 

The  results  will  be  presented  in  rectangular  cartesian  coordinates 
and  in  spherical  coordinates.  Some  simple  potential  functions  will  be 
analyzed. 

We  then  consider  methods  for  recovering  the  potential  function  and  the 

gravity  vector  from  the  gradient,  and  conclude  with  an  estimate  of  the 

necessary  sensitivity. 

* 

3.2.1  Gravity  Gradient  as  a  Function  of  the  Potential 

It  is  shown  in  classical  mechanics  that  the  forces  produced 
by  a  gravitating  body  can  be  represented  conveniently  by  means  of  operations 
on  a  scalar  point  function  of  position  called  the  potential  of  the  field.  If  the 
mass  distribution  is  given,  the  potential  of  the  field  can  be  obtained  by  adding 
the  separate  potentials  of  each  elementary  piece  of  mass  algebraically,  or 
equivalently,  by  integrating  over  the  volume.  Potentials  superimpose. 
Further  these  potentials  are  conservative.  By  this  is  meant  that  the  work 
required  to  move  a  mass  from  one  point  to  another  in  the  field  depends  only 
upon  the  initial  and  final  positions  and  not  upon  the  path  taken  between  the  two 
points.  For  a  unit  mass  this  work  is  numerically  equal  to  the  change  in 
potential  between  the  two  points.  If  two  points  have  the  same  potential  no  net 
work  is  required  to  move  from  one  point  to  the  other.  It  is  customary  to  take 
the  reference  level  of  potential  as  zero  at  infinity. 


10 


DS  64-R531-35 


The  potential  of  an  elementary  particle  of  mass  m  concen- 
I  trated  at  a  point  is 

Gm 
r 


V  =  - 


3.2-2 


r  is  the  distance  from  the  mass 
G  is  the  universal  gravitational  constant,  in  units 
appropriate  for  V,  m,  and  r. 


When  m  is  in  grams,  r  is  in  centimeters,  V  is  in  ergs,  and 
G  then  becomes  6.  67  X  10"  cm'Vgm/sec^.  This  number  is  not  known 
very  exactly,  although  the  associated  gravitational  acceleration  on  the 
earth  is  well  defined. 


We  calculate  the  gravity  force,  £ 
gradient  of  V.  In  traditional  vector  symbology 


by  taking  the  negative 


=  -  grad  V 

=  -  V  v. 


3.2-3 


In  this  notation  £  is  the  force  on  a  unit  mass  placed  in  the 
gravity  field.  For  the  particle  field  of  equation  (3.2-2) 


3.2-4 


We  shall  designate  unit  vectors  by  the  notation  1^  ,  meaning 
the  unit  vector  associated  with  the  q  coordinate.  Equation  3.2-5  is  just 
Newton's  equation  for  the  forces  of  attraction  between  a  unit  mass  and  a,  mass  m, 
inversely  proportional  to  the  square  of  their  separation  and  directed  toward 
each  other . 


For  more  general  potential  functions  we  may  explicitly  express 
£  as  a  function  of  the  partial  derivatives  of  V  with  respect  to  the  coordinates 
of  the  system  of  representation. 
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In  rectangular  cartesian  coordinates  we  have: 


*  lx.1*  *  *9*1 

In  polar  spherical  coordinates  we  have: 

'ft*  flfr'i****  rTfco 

•  if*fr  *'  1*  2 e  +  i‘4.  f  '4. 


3.2-6 


3.2-7 


3.2-8 


3.2-9 


Note  that  the  point  of  observation  is  not  moving  in  space. 

Thus  the  vector  shown  is  not  the  vector  which  would  be  indicated  by  a 
plumb  bob  on  the  earth.  It  is  the  gravity  vector  which  would  be  experienced 
by  either  a  satellite,  or  a  body  attached  to  a  non-rotating  planet. 

Following  our  definition  of  gravity  gradient  in  the  opening 
sentences  of  this  section,  we  can  evaluate  it  in  terms  of  partial  derivatives 
of  the  potential  function' and  the  direction  cosines  of  the  sensitive  axis. 

In  rectangular  coordinates  we  write  (gg)  in  matrix  form: 


(ffi  -  °T lv*;Kjl c 
=  C  ci  cz  cj) 


3.2rl0 


>*y.  a, 


ff  Ji£  JiiL  b2 

Ty* 

JlX.  J*X-  i*V.  C3 


3.2-11 


Note  that  |V*j£j  I  is  a  symmetric,  square  matrix.  The 
quantities  C,,  C2  an<*  C3  are  the  direction  cosines  of  the  sensitive  axis 
on  the  X,  Y,  ZT  axes  respectively.  Their  squares  add  to  unity. 


+  C22  +  c32 


3.2-12 
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Performing  the  indicated  operations  in  rectangular  coordinates 

pL  +  c»a  J&C 

°x  dy 


+  Cj 
2 


*yZ 


3.  2-13 


+  2C,CjpL  +2C,C3£V.  +2CzC3  IX 


x 


'by.'by 


In  spherical  coordinates  we  go  directly  to  the  extended  statement: 

*  .  f  Cjz  ces  G  Cr  \  J.  j 

/  w  r\  r  si 


r  sip  G 


cot  *)  r  Cm  <9  l^T 

+  C,1  -ift  +  -2C’-C«  +  *£r£i 

r  )V  *  ir  ^  r  sin  6 

df  *  r*  3/fl  s,n 


sa 


i$^r* 


3.2-14 


Before  continuing  we  want  to  note  some  other  relations.  An  invariant 
of  the  potential  field  is  the  divergence  of  the  gradient,  which  is  identically 
zero  wherever  there  is  no  mass.  This  is  Laplace's  equation. 


V^V  =  0 

In  rectangular  coordinates 


3.2-15 


if v 


alii  +ilv.  =  <> 


Hz 


3.2-16 


In  spherical  coordinates 


/ 


r2  sin  6 


(s,r*  &  & 


% 


£}■ 


3.2.  -17 

0 


If  the  potential  is  given  everywhere  on  a  closed  surface 
the  potential  function  is  uniquely  determined.  This  follows  from  the 
divergence  theorem  of  vector  calculus.  (Dirichlet's  problem). 

If  the  normal  derivative  of  the  potential  is  given  every¬ 
where  on  a  closed  surface  the  potential  function  is  uniquely  determined. 
(Neumann’s  problem). 
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Note  that  equations  (3.2-16,  ^7)  holding  in  the  mass -free 
space,  can  be  solved  by  classical  methods  to  provide  literal  solutions 
for  V,  and  correspondingly  for  its  gradient,  the  gravity  field.  Accordingly, 
if  such  a  solution  is  set  down,  and  if  the  literal  coefficients  can  be  evaluated 
from  measurements  of  (gg)  we  have  found  a  way  to  evaluate  £  outside  of  the 
body.  This  method  will  be  examined  below. 


Consider  a  gradiometer  at  sequential  points  on  a  trajectory, 
the  gradiometer  holding  a  fixed  attitude  in  X  Y  Z,  that  is,  inertial  space. 

If  the  gradiometer  is  oriented  along  the  X  axis,  is  unity  and  and  C 3 
are  zero.  Then 


3.2-18 


If  we  should  employ  a  fixed  orthogonal  triad  then  the  gradio- 
meters  would  measure  V_  and  a  check  for 

Tx*  >  tpr  >  fpr 

behavior  would  yield  a  sum  of  zero.  Note  that  the  reference  orientation 
is  arbitrary.  Therefore  any  inertially  fixed  triad  would  satisfy  this  con¬ 
dition,  regardless  of  the  trajectory,  or  the  coordinate  system. 


3.2.2  Analytic  Determination  of  the  Gravity  Field 


We  have  earlier  introduced  spherical  coordinates  because 
of  the  approximate  spherical  shape  of  planetary  objects,  and  the  equipotential 
surfaces  are  approximately  spherical.  Accordingly  we  expect  the  principal 
terms  in  an  analytical  description  of  the  field  to  be  particularly  simple  if 
we  use  this  coordinate  system. 

On  this  basis  a  literal  general  representation  of  the  potential 
field  in  the  space  outside  the  gravitating  body  can  be  obtained  as  the  sum  of 
terms  of  the  form 


V„  =  r 


_  -»-/ 


3. 2-19 


The  quantities  Yyj  are  functions  of  0  and  £  and  n  is  an  integer 
The  spherical  harmonics  Yu  can  be  written  as 


Yn  *  £  +  2.  Sm  m  4 

mmo  m  =0 


3.2-20 
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or,  using  complex  variables 

Yn 

m  M  -  n 


Thus  the  potential  field  is  given  by 

<»  n 


3.2-21 


V-  I  jV"''! 

7l*t  C  •»«-*  J 


3.2-22 


The  quantities  are  trigonometric  functions  of  0  alone,  and  are 

called  associated  Legendre  functions.  When  m  is  zero  they  are  called 
Legendre,  simple,  or  zonal  harmonics .  They  are  most  conveniently 
defined  by  introducing  a  new  variable  t. 


t  35  C .05  G 


3.2-23 


i  d"(tx~0 


^  =  -*V  dt  * 

EVf*-<y*  jmp* 


dt 


or 


3.2-24 


3.2-25 


3.2-26 


In  the  variable  t  these  functions  are  polynomials,  and 
t  ranges  over  the  interval  —  /  £.  t  /. 


Table  I,  taken  from  Heiskanen,  lists  some  of  the  lower  order  terms. 

The  special  term  is  called  the  sectorial  harmonic,  and  the  remaining 
terms  for  which  m  is  neither  zero  nor  n  are  called  tesseral  harmonics. 

We  shall  use  this  description  to  find  the  gravity  field. 
Suppose  we  have  determined  ^  ,  the  radial  gravity  gradient, 

at  several  positions  on  an  orbit.  r 
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Differentiating  our  equation  for  V  we  find: 

nr  =-fpr  =  3-2-27 

Thus  for  several  values  of  r,  Qt  $  we  have  values  for  ffr  ■  We 
number  each  data  point  to  distinguish  it.  Then  we  compute  the  value  of 

;>  w  y.  ~n-$  _  im  d 

jf  )  •  and  e  at  each  point  for  the  values  of  n  and  m 

we  consider  significant.  We  get  a  set  of  linear  equations. 

7fr: g  *  5"  ^mvz 

•  •  •  •  3.2-28 

•  •  •  • 

•  • 

a  a  ,  —  #  # 

In  these  equations  the  k  &&}  are  known,  and  the  ffy  terms  are  known. 
The  0,^  terms  are  not  known.  Accordingly,  we  can  solve  the  equations 
(invert  the  matrix)  if  we  have  a  sufficient  number  of  data  points. 
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Further,  once  we  have  found  the  O rtm  y  we  can  find  the  radial  com¬ 
ponent  of  the  gravity  vector  by  the  relation 

gr  *  -  +  ^  r"n  2^Cnm  pnme  1  m  ^  3.2-29 

*r 

The  other  two  components  are  similarly  available. 

If  we  have  more  data  than  necessary,  for  the  number  of  terms  considered 
significant,  we  can  make  use  of  the  redundancy  to  provide  a  least  squares 
fit  to  the  data. 

By  this  means  we  can  obtain  a  global  description  of  the 
potential  function  and  the  gravity  vector  everywhere  outside  the  planet  from 
which  we  can  get  a  reference  surface  or  selenoid  or  geoid. 

3.2.3  A  Method  for  Mapping  Anomalies 

Potential  functions,  at  least  of  gravity  fields,  are  linear 
field  properties.  This  is  a  characteristic  of  space.  Accordingly,  if  we 
have  a  potential  field  we  may  partition  it  into  several  potential  fields,  since 
the  sum  of  potential  fields  is  itself  a  potential. 

This  property  is  useful,  for  example,  in  handling  the  reduction 
of  the  gradient  data  into  gravity  anomalies  on  a  reference  geoid.  It  has 
appeared  implicity  in  the  data  reduction  process  described  above,  where 
we  have  taken  a  sum  of  functions  (the  spherical  harmonics  multiplied  by 
appropriate  powers  of  the  radius)  to  represent  the  total  potential. 

Each  of  the  terms  T  *  *  vw  is  a  potential  function  in  its  own 
right. 


Assume  that  a  satisfactory  geoid  has  been  found,  for  example, 
by  following  the  method  of  the  previous  section.  For  this  discussion  we 
assume  further  that  it  accounts  for  the  first  few  terms  of  the  total  potential. 

If  we  compute  the  gravity  gradient  associated  with  this 
reference  potential  and  subtract  it  from  the  data,  the  residuals  constitute 
the  anomalous  gravity  gradient  data. 

A  characteristic  of  the  series  representation  we  have  used 
is  that  the  fine  (anomalous)  structure  of  the  field  requires  many  terms  in 
the  expansion  to  provide  a  satisfactory  approximation.  However,  we  do  not 
really  want  the  analytic  description  for  most  of  the  interpretive  work  of 
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anomaly  identification.  What  is  required  is  a  map  of  the  anomalous 
behavior  on  the  reference  surface,  in  a  local  manner. 

In  the  following  paragraphs  we  will  set  down  a  theory. 

We  begin  by  returning  to  fixed  rectangular  coordinates 
and  the  properties  of  the  potential  function  in  this  system. 


Let  V  be  a  scalar  potential  function.  Then  its  partial 

jv  ,2X-  ,2X- 

derivatives  Jx  which  are  the  components  of  the  gravity 

field  vector  are  themselves  potential  functions.  They  satisfy  Laplace's 
equation.  This  is  demonstrated  readily  by  direct  computation. 


JVvJ 


+32L 

f 


V 

[oj 


3.2-30 


...  3  2-31 

We  have  used  the  commutative  property  of  partial  differ¬ 
entiation  to  arrive  at  the  result. 


Now,  suppose  that  we  had  measured  the  derivative  of  gx 
in  the  x  direction  (the  gravity  gradient)  on  some  plane  on  which  x  is  constant 
That  is,  we  know  as  a  function  of  y  and  z  everywhere  on  the  plane. 

We  can  now  use  another  relation  of  potential  theory  to  find  gx  itself  on  some 
other  plane.  The  relation  assumes  we  are  concerned  with  two  potential 
functions  U  and  W,  in  mass-free  space.  Then,  On  a  closed  surface  S, 


\  (u¥t  ~wj^)As  3 

Think  of  W  as  the  $x.  field.  Then  zv/dx  is  the 


3.  2-32 


measured  gravity  gradient  if  part  of  the  closed  surface  is  the  x  plane.  On 
this  plane  we  would  like  )tf/ bK  (the  normal  derivative  of  U)  to  be  zero. 
We  accomplish  this  by  an  educated  guess.  We  take  U  as  the  function 
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U  *  [C«  +  (f-y)1*  &  ~P  2J  3.2-33 

+  Jjf«f  3 c/  *  (/9-fi2*  &  "&  j 

U  it  now  a  function  of  6  quantities  ;  7“  and  x,  y,  z. 

At  the  point  x  =  Ot  'j  - yS  z  -  ^  in  the  volume  above 

the  x  =  0  plane,  U  is  singular,  and  we  must  somehow  exclude  it  from  the 
•pace  surrounded  by  S.  We  do  this  by  imagining  a  small  sphere  S'  drawn 
around  this  point  and  attach  this  surface  to  S;  the  whole  surface  over 
which  we  integrate  becomes  S  +  S'..  The  normal  to  S  is  in  the  negative  x 
direction,  the  normal  to  the  sphere  is  into  the  sphere. 

d</  / 

Since  on  the  x  plane  /  3x  is  zero,  we  have 

5  5  (•'£-*&)*' = ° 

3  f  9^,/  /  / 

Over  the  small  sphere,  U  is  constant  in  value  and  i _ e<S 

0  uM  W  ^ 

is  zero.  On  the  other  hand,  0^,  w  over  the  sphere  S'  is  essentially 

equal  to  the  constant  value  of  W  at  the  point  (  <*,  }  )  and  ^  ^  • 

is  4  V  .  ■ 


3.2-34 


(  d\J 


Therefore,  the  final  equation  is 


4 ir  WC*-,/?,  r)  =  -  ^ 


3.2-35 


♦*-?x 


Now  taking  for  W  we  have: 

o  f  [(^'/  +  C/f-p2  +  * 


3.2-36 


Note  that  if,  by  definition,  we  take  the  reference  plane  as 
x  =  0  the  result  becomes  more  compact. 


Sxm~  J  [°‘2+(/?-y)2+(i~£)zJ 


3.2-37 
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Practically,  we  do  not  know  $$  x  continuously  over  the 
plane.  At  best  we  know  it  in  terms  of  samples;  and  a  discrete  numerical 
integration  process  is  taken  to  approximate  the  above  operations. 

This  process  is  exactly  that  evaluated  by  Paterson  in 
Geophysics,  August  *61.  In  the  article  Paterson  deals  with  the  problem 
of  converting  an  airborne  survey  of  gravity  gradient  onto  the  nearby  earth 
as  gravity  force,  and  proposes  a  numerical  technique  for  making  the 
computations . 

This  appears  to  be  a  reasonable  approach  to  mapping  the 
anomaly  field.  However  further  work  is  needed,  since  the  expected  satellite 
observations  will  not  be  on  a  plane  surface.  It  is  suggested  that  equivalent 
theory  for  a  spherical  coordinate  system  be  developed  along  with  appro¬ 
priate  numerical  methods.  As  an  alternative  the  spherical  harmonic  repre¬ 
sentation  can  always  be  used. 

3.2.4  An  Estimate  of  Required  Instrument  Sensitivity 

We  now  have  a  data  reduction  procedure  from  which  we  can 
develop  a  reference  potential  surface  (geoid)  and  an  approach  to  locally 
mapping  the  residual  gravity  field,  the  anomaly  field,  onto  the  geoid. 

In  order  further  to  evaluate  the  instrumental  problem,  we 
must  develop  some  estimates  of  the  necessary  sensitivity  of  the  gradiometer. 
We  would  like  to  relate  this  to  the  more  familiar  gravity  anomaly  observed 
on  the  planetary  surface.  To  do  this  we  have  employed  a  simple  anomaly 
model,  consisting  of  a  large  central  planetary  mass  M,  and  a  small  mass  m 
Ibcated  near  the  surface  of  the  planet.  For  an  orbit  at  100  miles  above  the 
planet  we  have  calculated  the  radial  gravity  gradient  above  the  anomaly  and 
plotted  this  against  the  surface  gravity  anomaly,  for  various  sites  of  the 
anomaly  below  the  surface. 

For  this  model  we  have 

V  =  -  {  GM  +  Gm  )  3.2-38 

R  R-rc  ' 

M  =  planetary  mass;  for  the  earth  around  10^^ 
tons;  see  table  II  for  other  planets. 
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xn  =  excess  mass  associated  with  the  anomaly. 

For  visualization  note  that  one  cubic  mile  at 
mean  earth  density  weighs  10^  tons. 

R  =  altitude  of  orbit  above  planet  center. 

r0  =  anomaly  site  above  planetary  center. 

V  =  potential  along  line  joining  the  two  masses, 
beyond  the  surface. 

h  =  R  -  rD;  altitude  of  orbit  above  the  anomaly. 
hQ  =  altitude  of  surface  above  the  anomaly. 


Radial  gravity  =  gr  = 


Radial  gravity  gradient  =  gg  =  2  -pM_  +  ..iP.™- 

R3  (R-rc)3 


3.2- 39 

3.2- 40 

3.2- 41 

3. 2- 42 


In  these  equations  if  R  is  very  closely  the  planetary  radius 
then  GM/R4  is  the  normal  surface  gravity  and  2GM/R^  is  the  normal 
surface  gravity  gradient.  Further  m  f  R  \  2  is  then  the  fraction  of 

M  \h/ 

normal  gravity  due  to  the  anomaly  and  m  /  RV  is  the  fraction  of  normal 

M  \  hj 

gravity  gradient  due  to  the  anomaly.  This  latter  defines  the  necessary 
effective  sensitivity  of  the  instrument.  In  the  sense  used  here  sensitivity 
is  more  often  called  "dynamic  range",  the  ratio  of  the  least  detectible 
measurement  to  the  measurement  range.  We  define  as  the  discrimination 
level  the  value  of  gravity  gradient  corresponding  to  S,  Thus  if  the  range  of 
measurement  is  10'^  g/foot  and  S  is  10  the  discrimination  level  is  10-^2g/foot. 
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In  accord  with  this  philosophy  we  now  can  derive  a 
relation  between  the  surface  gravity  anomaly  and  the  corresponding 
change  in  gravity  gradient  at  altitude.  At  the  surface,  the  ratio 
anomalous  to  normal  gravity  is: 


3.2-43 


At  altitude  the  ratio  of  anomalous  to  normal  gradient  is^ 


)r 

-JO-  / 

r*-' 

’(¥)W 

)r 

M  ' 

\  h  / 

'  M 

1  he' 

\  J  \  h  / 

3.2-44 


h 


r-o 


3.2-45 


h3 

On  earth  gr  is  1  gravity  a:  10^  milligals.  For  an 


anomalous  body  or  mass  one  mile  below  surface,  producing  a  one  milligal 
gravity  deviation  at  the  surface  and  for  a  100  mile  altitude  orbit  S  is 
4  x  10-9 


S  has  been  calculated  for  earth  and  lunar  orbits,  taking 
lunar  gravity  as  160,  000  milligals,  and  lunar  radius  as  1000  miles.  The 
important  results  are  tabulated  below. 


Earth  Moon 


milligals 

ho 

=  1  mile 

5  miles 

1 

mile 

5 

miles 

1 

4  x 

10'9 

10“7 

6 

X 

10-8 

1  .  6 

x  10'6 

10 

4  x 

10-8 

10‘6 

6 

X 

io-7 

1.5 

x  10'5 

100 

4  x 

10'7 

10'5 

6 

X 

10‘6 

1.5 

x  10"4 

1000 

4  x 

10"6 

10'4 

6 

X 

IO"5 

1.5 

x  10"3 

Value  s  of  S 
Table  III 
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It  appears  that  for  a  system  with  a  sensitivity  S  of  10"^  the 
system  would  indicate  lunar  anomalies  of  one  to  10  milligals  at  the  surface, 
arising  from  anomalous  structure  one  to  five  miles  below  the  nominal 
planetary  surface.  At  a  level  of  10"®  the  surface  anomalies  would 
correspond  to  10  to  100  milligals.  This  latter  is  considered  to  be  a  rea¬ 
sonable  level  to  start  the  investigation  of  the  lunar  gravity  field  from  an 
orbiting  satellite.  The  discrimination  level  for  an  S  of  10"^,  in  a  I0“7g/foot 
field  is  simply  10"^g/foot,  or  0.03  E  (E  otvos  units).  There  is  no  question 
but  that  this  is  an  extremely  delicate  measurement,  and  potentially  vul¬ 
nerable  to  a  vast  variety  of  normally  hidden  noises.  For  example  the  potential 
field  of  a  one  kilogram  mass  of  mean  earth  density  material,  at  a  distance 
of  one  foot  produces  a  gravity  gradient  500  times  as  great  as  that  we  seek 
to  measure.  A  shifting  of  this  mass  by  0.  1  mm  can  result  in  a  signal  change 
approximating  the  discrimination  level.  But  this  can  be  interpreted  as  a 
constraint  on  laboratory  design.  For  example  the  gradiometer  could  be 
put  on  a  boom  so  it  is  remote  from  the  major  mass.  Then,  mass  motions 
within  the  satellite  will  not  affect  the  gradiometer  output. 

The  instrument  must  also  be  protected  from  developing 
thermal  gradients  along  its  length  since  this  will  result  in  a  non- symmetrical 
configuration.  If  the  end  string  lengths  differ  due  to  unequal  expansion  a 
biased  output  will  result.  Such  a  disturbance  must  be  kept  below  the  dis¬ 
crimination  level. 

Periodic,  automatic  calibrations  of  the  instrument  will  allow 
for  the  compensation  of  temporal  shifts  in  the  instrument  scale  and  bias. 

These  shifts  are  unavoidable  in  any  delicate  instrument  over  a  long  period 
of  time.  It  is  reasonable  to  expect,  however,  that  such  shifts  can  be  de¬ 
tected  to  the  discrimination  level  and  compensated  systematically. 

It  is  also  important  to  note  that  the  discrimination  level  is 
reduced  by  virtue  of  the  repetitious  data  collection  inherent  in  the  orbiting 
measuring  system.  Random  noise  will  be  averaged  to  a  low  level  by 
statistical  techniques  included  in  the  data  processing. 

These  problems  of  resolution,  noise,  thermal  effects  and 
statistical  data  processing  have  already  been  solved  for  the  Arma  Vibrating 
String  Accelerometer. 

Considering  the  foregoing,  the  dynamic  range  of  the  vibrating 
string  gradiometer  can  be  set  so  that  the  discrimination  level,  correspond¬ 
ing  to  lunar  surface  anomalies  of  10  -  100  milligals  will  be  achieved  while 
the  maximum  gradiometer  output  will  be  at  the  level  of  the  normal  lunar 
gradient. 


1 
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3.2.5  Effects  of  the  Sun  and  the  Earth 

Two  possibly  important  planetary  ambient  sources  should 
be  noted;  those  due  to  the  earth  and  the  sun.  In  computation  of  the  lunar 
field  these  can  and  should  be  subtracted  from  the  gradient  data  before 
processing. 


The  gradient  field  of  the  earth  in  the  vicinity  of  the  moon 
is  approximately  4  x  10"  g/foot  and  varies  about  50%,  from  about 
8  x  10  to  4  x  10"  13. 

The  gradient  field  of  the  sun  is  approximately  2  x  10“ 
g/foot  in  the  vicinity  of  the  moon.  This  appears  ignorable  because  it  is 
two  magnitudes  less  than  the  discrimination  level  we  have  been  discussing. 

Note  that  the  attitude  of  the  gradiometer  relative  to  the 
radius  to  the  earth  must  be  entered  into  the  corrective  computation.  This 
can  be  determined  readily  from  the  attitude  of  the  sensor  with  respect  to 
the  moon. 


3.2.6  Summary 

In  this  section  we  have  presented  the  basic  gravity  gradient 
theory  and  a  method  for  converting  gradient  observations  into  gravity  field. 

An  analysis  of  the  anomaly  field  of  a  point  anomaly  indicates  that  a  system 
sensitivity. of  one  part  in  105  is  desirable  for  anomaly  mapping.  It  is 
further  noted  that  at  this  level  account  must  be  taken  of  the  gradient  field 
of  the  earth.  This  is  a  simple  computation. 

3,  3  Indications  of  a  Gradiometer  on  a  Rotating  Platform 

A  useful  physical  picture  of  the  gradient  phenomenon  is  obtained  if 
one  imagines  two  masses  fastened  together  by  an  ideal,  short,  infinitely 
slender  rod  and  placed  in  the  field  of  a  gravitating  body.  If  the  rod  is  oriented 
radially,  the  gravitational  force  on  the  inner  mass  is  greater  than  that  on 
the  outer  and  the  rod  is  in  tension  by  an  amount  equal  to  the  radial  gradient 
of  the  gravity  force  field  multiplied  by  the  length  of  the  rod.  If  the  rod  is 
now  held  horizontally,  the  convergence  of  the  acceleration  forces  toward 
the  center  of  the  gravitating  body  results  in  a  force  along  the  rod  tending 
to  bring  the  masses  together,  a  compression.  It  is  clear  that  this  is  a 
measure  of  the  curvature  of  the  gravity  field  at  the  point.  Numerically,  for 
the  simple  force  field  assumed,  the  compression  is  one  half  the  tension. 

As  the  rod  is  turned,  it  passes  from  tension  to  compression  and  back  again 
in  one  half  turn. 
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Now  visualise  the  same  experiment  performed  on  an  orbiting 
satellite.  To  be  specific,  assume  a  circular  orbit.  If  the  rod  -  the 
gradiometer  -  is  held  fixed  in  inertial  space,  its  indications  will  be 
independent  of  the  orientation  of  the  satellite,  but  will  vary  with  time 
fcince  the  orientation  of  the  rod  with  respect  to  the  radius  to  the  planet 
varies  cyclically  in  orbit.  Thus,  the  rod  will  sample  the  gravity  gradient 
in  various  directions,  as  it  did  in  the  non-orbiting  experiement  described 
above.  If  the  gradiometer  should  be  continuously  controlled  to  lie  along 
the  radius,  then  it  will  measure  the  radial  gradient  plus  the  square  of  the 
orbital  angular  velocity  in  inertial  space.  Because  of  the  simple  orbital 
geometry  this  angular  term  is  exactly  one  half  the  gradient  term  and  has 
the  sense  of  tension.  A  gradiometer  in  the  orbital  plane  lying  normal  to 
the  radius  would  read  zero  (its  compression  is  reduced  by  the  equal 
tension  due  to  the  angular  velocity).  A  gradiometer  perpendicular  to  the 
orbital  plane  would  read  the  usual  compressive  force  due  to  gravity  alone. 
Accordingly  with  the  addition  of  rate  gyroscopes,  the  three  components 
of  gravity  gradient  can  be  obtained  separately.  It  can  be  shown  that  the 
sum  of  the  three  quantities  measured  by  a  triad  of  gradiometers  is  twice 
the  square  of  the  angular  velocity  of  the  gradiometers  in  inertial  space. 


If  the  gradiometer  is  on  a  rotating  platform  its  indications  are 
composed  of  two  parts;  the  normal  gravity  gradient  existing  at  the  point 
in  space,  in  the  direction  of  the  sensing  axis  plus  a  component  due  to 
angular  velocity  of  the  platform.  If  we  follow  the  definition  of  gradient 
employed  earlier  we  find  that  the  term  due  to  angular  velocity  is  Wg2  where 
JWgJ  is  the  magnitude  of  the  component  of  angular  velocity  normal  to  the 
sensing  axis.  We  compute  Wg  from  the  vector  relation: 


W  =  1  x  W 

8  8  - 


3.3-1 


The  sum  of  the  squares  of  the  three  Wg  terms  is  equal  to  2W2. 
Accordingly,  the  sum  of  the  3  gradiometer  outputs  on  a  rotating  platform 
corresponds  to: 

88j  "*■  882  +  883  =  2W^  -  .  In  mass-free  space  is  zero. 

If  a  gradiometer  orthogonal  triad  is  constrained  so  that  one  gradio¬ 
meter  tracks  the  geocentric  radius  and  one  tracks  the  local  horizontal 
along  the  trajectory  then  the  angular  velocity  of  the  platform  is  determined 
by  the  orbital  parameters.  For  a  simple  central  force  field  and  a  circular 
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orbit  the  results  are  especially  simple  and  useful  for  making  system 
estimates.  Denoting  indications  by  I  we  have: 


if  =  I  ffr 


(radial)  3.3-2 


Jg  —  0  (along  trajectory)  3.3-3 

and  ^r-//r=-2^2  3.3-5 

The  square  of  the  orbital  angular  velocity  is  one  half  of  the  radial 
gravity  gradient. 


(normal  to  3.3-4 

orbit  plane) 


3.4  Instrumentation  for  Gradiometry 

Several  instruments  which  may  serve  as  gradiometers  are  described 
below.  Each  instrument  measures  the  difference  between  the  forces 
affecting  two  (or  more)  proof  masses  when  these  masses  are  separated 
along  the  sensitive  axis.  The  different  proof  masses  are  self  evident  in  the 
discussion  of  matched  accelerometers,  the  vibrating  string  gradiometer 
and  the  torsional  gradiometer.  For  the  accelerometer-on-a-boom  and  the 
San  Marco  experiment  the  orbiting  satellite  itself  is  one  of  the  proof  masses. 
The  gas  diffusion  equipment  is  made  up  of  innumerable  proof  masses  whose 
relative  motion  may  be  analyzed  statistically. 

The  gradiometer  which  appears  most  suitable  for  the  above  described 
measurements,  is  an  extension  of  equipment  already  developed  at  the  Arma 
Division.  This  is  the  vibrating  string  accelerometer  which  has  been  used 
so  successfully  in  inertial  guidance  systems  and  which  is  also  the  basis  for 
ruggedized  thrust  instrumentation  recently  delivered  to  NASA- Lewis  Research 
Center.  A  low  range  accelerometer  for  ion  propulsion  measurements  is  now 
under  evaluation  for  RTD. 
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3.4.1  Basic  Theory  of  the  Vibrating  String  Instrumentation 
Single  String,  Single  Mass  System 

A  string  under  tension  will  vibrate  when  excited.  The  fre¬ 
quency  of  vibration  is  a  function  of  the  dimensions  of  the  string,  the 
string  density  and  the  tension  of  the  string.  The  frequency  of  vibration  (f) 
is  directly  proportional  to  the  square  root  of  the  tension  of  the  string  (T) 
and  inversely  proportional  to  the  square  root  of  the  mass  of  the  string  (m) 
and  its  length  (L). 

f  =  1/2  T_  3. 4-1 

\  mL 

For  a  string  supporting  a  mass  M  the  tension  on  the  string  is  equal  to  the 
initial  tension  (T0),  due  to  the  mass  attraction  of  the  gravitational  field, 
plus  any  other  force  acting  on  the  sensitive  mass. 

T-T0  +  Ma  3.4-2 

In  this  case  TQ  is  Mg  and  the  inertial  reaction  force  is  Ma 


UlLLLLU 


String 


Proof  Mass 


Mounting 

Frame 


From  this  it  can  be  seen  that  for  any^particular  string  and  mass,  the 
frequency  of  vibration  can  be  varied  by  changing  the  acceleration  applied 
to  the  sensitive  mass.  If  we  consider  the  acceleration  as  variable,  the 
frequency-tension  relationship  can  be  expressed  as  a  Taylor  Series  ex¬ 
pansion  about  zero  acceleration. 


or 


P=  K0  +  K,*+  k.o.  3-  +  K3o-3  + 


3.  4-3 
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When  a,  the  applied  acceleration,  is  zero,  f  =  Kq  which  is  called  the 
bias  frequency.  With  applied  acceleration,  there  is  a  linear  change  in 
frequency  plus  other  non-linear  terms. 


Double  String,  Single  Mass  Systems 


The  even  order  non-linear  terms  in  the  Taylor  Series 
expansion  can  be  reduced  in  magnitude  by  using  two  strings  along  one 
axis  with  the  mass  between  them. 


Mounting 

Frame 


In  this  representation,  the  coefficients  approach  zero  as  tl  '  c'  der  of 
the  term  approaches  infinity  and  the  series  converges.  If  we  express 
the  frequency  of  each  string  as  a  series,  we  have: 


f,  *  Ko,  +K„°-  +  Kii«-2  +  K3i°-,+  --.  For  T=T„  +M0./2 
$2  =  ^03.  ~  K I2a'  *^22al  ~ ^12a',+  '  *  ■  1~=T0  -  Mo. /2. 

3.4-4 


Since  the  K.  terms  are  dominant,  the  frequency  of  one  string  will  increase 
while  the  otner  will  decrease.  Taking  the  difference  of  these  two  fre¬ 
quencies  we  have: 

frFx=(K 0,  -  x;a  )*<X + «a)  «  +  04/  -  +  Kx) «3+  ■■■ 


3.4-5 


If  the  two  strings  are  closely  matched,  all  the  even  terms  in  the  difference 
frequency  equation  become  negligibly  small  and  only  the  odd  terms  remain. 
Since  the  series  converges  rapidly,  the  only  non-linear  term  of  any  signi¬ 
ficance  is  the  K3  term  but  this  can  be  compensated  if  a  is  known  approxi¬ 
mately.  This  method  is  useful  in  achieving  a  substantially  linear  device. 
However,  because  the  two  strings  are  attached  to  a  single  mass,  there 
may  be  coupling  between  the  strings  when  they  are  near  the  same  frequency. 
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The  ameliorative  i«  shown  below. 

Doable  String,  Double  Mass  System 


Two  modifications  are  made  on  the  double  string,  single 
mass  system  to  minimize  the  effect  of  coupling  between  the  strings. 
Two  separate  but  equal  sensing  masses,  joined  with  a  relatively  soft 
spring,  are  used  in  place  of  the  single  mass.  ^ 

Strthf  ____ 


#/ 


kTJWSh 


Strt*} 

*2 

• 

Mass  *2  \ 

'  Frame 


^  Mass  #/ 

This  reduces  the  transmission  of  energy  from  one  string  to  the  other  while 
not  causing  any  deterioration  of  the  system.  The  second  modification  is  the 
substitution  of  flat  tapes  for  the  strings:  These  tapes  are  oriented  to 
vibrate  in  mutaully  orthogonal  planes  and  thus  reduce  the  tendency  to 
vibrate  sympathetically. 


Cross  Centering  System 

In  order  to  reduce  the  effect  of  accelerations  not  directed 
along  the  axis  of  the  tapes,  it  becomes  necessary  to  restrain  the  masses. 
This  is  accomplished  by  cross  supports  placed  perpendicular  to  the 
sensitive  axis.  The  following  diagram  shows  the  basic  Arma  Vibrating 
String  Accelerometer. 


Mounting  Flange 


’iWs 9S 

■ 

k 

\WHKM 

1 

Electrical  Considerations 


Permanent  magnets  provide  a  means  for  laterally  vibrating 
the  tapes  at  their  natural  frequencies.  If  a  current  is  passed  through  the 
tape  in  a  magnetic  field,  a  force  is  produced.  This  force ^tends  to  move 
the  tape.  If  the  current  is  reversed,  the  force  is  in  the  opposite  direction 
and  the  tape  is  moved  in  that  direction. 
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In  this  manner,  the  tapes  can  be  vibrated.  AC  voltages  are  induced  in 
the  tapes  by  their  motion  in  the  magnetic  fields.  These  voltages  are 
regenerated  through  stable  high  gain  amplifiers  and  returned  to  the  tapes 
to  provide  the  energy  necessary  to  sustain  oscillation.  The  tape  is  there¬ 
fore  part  of  a  tank  circuit  which  oscillates  at  its  natural  frequency. 


The  dynamics  of  the  vibrating  tapes  in  the  magnetic  fields  provided  by 
the  pefmanent  magnets  have  properties  which  are  completely  analogous 
to  those  of  a  parallel  resonant  electrical  circuit.  This  feature  permits 
the  design  of  an  appropriately  controlled  feedback  amplifier  to  maintain 
a  constant  amplitude  of  oscillation. 

v 

Accelerometer  Operation 

The  two  proof  masses  are  free  to  move  along  the  sensitive 
axis  and  are  restrained  in  motion  along  the  cross  sensitive  axes  by  the 
suspension  system.  An  acceleration  along  the  axis  of  the  end  tapes  causes 
the  masses  to  increase  the  tension  of  one  tape  and  decrease  the  tension  of 
the  other  tape  by  approximately  an  equal  amount.  Since  the  frequency  of 
tape  oscillation  is  proportional  to  the  square  root  of  tension,  the  instan¬ 
taneous  frequency  difference  of  the  tape  oscillations  is  a  measure  of  the 
acceleration. 


Because  of  the  square  root  relationship,  the  frequency 
difference  is  only  approximately  proportional  to  acceleration.  The  small 
effects  introduced  by  higher  order  terms  in  an  actual  accelerometer  can 
be  compensated. 

Vibrating  String  Gradiometer 

It  is  a  simple  step  from  the  Vibrating  String  Accelerometer 
to  the  Vibrating  String  Gradiometer.  A  third  Vibrating  tape  is  inserted 
between  the  two  proof  masses  and  its  tension  variation  will  be  a  function 
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of  the  difference  in  the  forces  affecting  the  two  masses,  in  particular 
the  gravity  gradient.  A  simple  sketch  of  the  gradiometer  is  shown 
below. 


L.  =  string  length 
M  =  proof  mass 
p  -  pole  magnets 
a  =  applied  acceleration 


e  =  incremental  acceleration 
=  L<2  x  gravity  gradient 


Simplified  sketch  of  Proposed  Gradiometer 


The  gradiometer  consists  of  two  proof  masses  M  held 
together  along  the  sensitive  axis  by  a  long  metal  tape  (string)  and  this 
combination  held  to  the  frame  by  two  short  end  strings  Lpand  L-j.  Each 
string  is  caused  to  vibrate  by  passing  an  alternating  current  through  it 
as  it  lies  in  a  permanent  magnetic  field.  The  frequency  of  that  current  is 
equal  to  the  natural  mechanical  frequency  of  the  string  and  is  established 
and  maintained  electrically  through  a  feedback  loop. 

The  natural  frequency  of  the  string  is  a  function  of  its  length, 
mass  and  the  tension  to  which  it  is  subjected.  When  the  unit  is  held  parallel 
to  the  local  vertical,  the  lower  proof  mass  is  subjected  to  a  greater  gravi¬ 
tational  force  than  the  upper  proof  mass.  Hence,  the  lower  end  string 
tension  is  reduced,  the  upper  end  string  is  increased,;  while  the  center 
string  is  increased  by  the  difference,  e,  in  the  gravitational  forces  on  the 
two  proof  masses.  This  difference  (gradient  times  length)  is  obtained  by 
processing  the  frequency  data  from  the  three  strings.  A  brief  mathematical 
description  of  the  relationship  between  the  instrument  parameters,  the 
differential  (gradient)  acceleration,  e,  and  the  ambient  acceleration,  a, 
follows . 
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where  f  =  natural  mechanical  frequency 

L  =  length  between  clamps 

T  =  tension  (force) 

=  linear  mass  density  (mass  per  unit  length) 

We  will  apply  this  to  the  two  mass,  three  string  system 
shown  in  the  figure  below.  Slightly  different  accelerations  are  applied  to 
each  proof  mass.  The  tensions  in  each  string,  which  were  each  T0  initially, 
change  in  such  a  manner  as  to  balance  the  externally  applied  forces  without 
changing  the  overall  length  of  the  system,. 


Therefore 


T 1  =  M(a-e/2)+T2 

3.  4-7 

T2  =  M(a+e  /2)+T  3 

3.  4-8 

Li  +  L2  +  L3  Lt  ;  A  Li  =  0 

3.  4-9 

These  changes  in  tension  are  the  new  values  of  tension  less  T0. 
in  length  are  directly  related  to  the  compliance,  S,  of  the  strings 
change  and  the  original  length. 

The  changes 
,  the  tension 

A  Lj  =  S  ^  Ti  Li 

3.4-10 

Then 

^  Ti  =  2Ma  +  AT3 

3.4-11 

ATz  =  Ma  +  Me/2  +  <^T3 

3.  4-12 
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from  equations  (3.  4-7,  8)  and 

Ll  A  Ti  +  L2  A  T2  +  L3  A  T3  =  0  3.4-13 

from  equation  (3.4-10) 

Solving  these  last  three  equations  for  the  string  tension  changes  we  find: 

ATl  =  -  hi  _Ma  3.4-14 

LT  T” 

A  T2  =  Me  3.4-15 

Lt 

A  T3  =  -ii.  £*«.  +Ma  3.4-16 


The  primary  quantity  of  interest  is  e,  the  gravity  gradient  multiplied  by  the 
length  of  the  second  string,  L.2.  Note  that  a  is  the  measure  of  the  accelera¬ 
tion  at  the  center  of  the  instrument.  If  the  instrument  were  on  a  satellite  in 
orbit,  a  would  be  the  gravity  gradient  at  orbit  altitude  multiplied  by  the  dis¬ 
tance  between  the  satellite  center  of  mass  and  the  instrument.  The  differ¬ 
ence  in  end  string  tensions  give  the  ambient  acceleration  while  the  tension  of 
the  second  string  gives  the  gravity  gradient. 


The  string  frequencies  are  related  to  the  instrument  parameters  and  acceler¬ 
ation  environment  through  the  following  equations: 


(r*-tT 


3.  4-17 


3.  4-18 


3.  4-19 


These  assume  that  the  end  stringB  are  of  equal  length  and  that / O  ,  the  linear 
mass  density  is  the  same  for  each  string.  Further,  through  the  choice  of 
string  material  and  configuration  the  percentage  change  in  length  is  made 
negligible  compared  to  the  percentage  change  in  tension  from  eir  respective 
nominal  values.  Then  the  frequency  depends  only  on  the  change  in  tension 
and  not  on  the  associated  change  in  length. 

If  the  frequency  of  the  middle  string  is  defined  as  fQ>  when  e  is  zero,  the 
gravity  gradient  is  given  by: 
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.  i  -ff i-  f  Jf.  LsJlL 

(}9)  -  LX  -  V*  fo  S  *  Ll 


The  ambient  acceleration  is  given  by: 


3. 4-20 


3.  4-21 


In  order  to  extract  the  gravity  gradient  from  equation  (3.  4-20)  it  is  necessary 
to  know  the  nominal  (bias)  frequency  f0  corresponding  to  the  initial  tension  Tc. 
This  frequency  will  be  set  when  the  instrument  is  manufactured  and  will  be 
monitored  as  part  of  the  periodic  calibration  procedure.  It  is  possible  also 
to  monitor  this  value  as  an  indication  of  instrument  stability  during  its  normal 
operation. 


For  this  data  processing  function  consider  the  relation: 


3.  4-21 


This  equation  provides  a  measure  of  T0,  and  can  be  used  to  test  for  off- 
nominal  behavior. 

In  summary,  it  has  been  shown  that  by  appropriate  combinations  of  the  three 
output  frequencies,  the  gravity  gradient  and  ambient  acceleration  can  be 
obtained  simultaneously.  In  addition,  a  monitoring  technique  is  available  to 
indicate  whether  or  not  the  instrument  is  in  need  of  re-calibration. 


If  the  environment  requires  it,  a  thermal  control  may  be  added  to  the  gradio- 
meter.  Also,  with  regard  to  operating  conditions,  if  it  is  necessary  to  store 
the  data  for  future  processing,  it  should  be  noted  that  the  frequency  modulated 
nature  of  the  output  signals  lend  themselves  readily  to  tape  recording  with¬ 
out  intermediate  coding  steps. 
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3.  4.  2  Accelerometers  for  Gradiomctry 

A  pair  of  high  quality  accelerometers,  separated  by 
an  accurately  known  distance  and  having  their  outputs  in  opposition,  give  the 
measure  of  the  gravity  gradient  over  that  distance  (taking  orientation  relative 
to  the  radius  vector  into  account).  It  may  be  possible,  from  an  orbiting 
satellite,  to  obtain  separations  between  the  accelerometers  of  the  order  of 
hundreds  of  feet. 

Suppose  a  100  foot  separation.  Then  a  10"5g  accel¬ 
erometer  could  be  used  to  measure  10“?g/foot  of  gravity  gradient.  The 
accelerometer  range  is  obviously  shifted  into  a  more  desirable  region.  How¬ 
ever,  this  does  not  improve  the  dynamic  range,  S,  of  the  instrument,  since 
the  threshold  is  generally  proportional  to  the  maximum  reading  for  a  given 
instrument  design  and  environment. 

In  addition  the  designer  is  now  faced  with  the  problem 
of  matching  performance  and  stability  of  two  instruments  over  their  entire 
operating  range,  in  the  face  of  different  thermal  environments.  The  separ¬ 
ation  must  be  continuously  measured  to  high  accuracy,  and  the  operational 
problems  of  erecting  and  holding  a  long  boom  in  different  attitudes  is  un- 
warrantedly  difficult. 

3.  4.  3  Gradiometer  on  a  Boom 


One  can  obtain  the  geometric  measure  required  for  the 
field  strength  and  gravitational  acceleration  determination  without  referring 
to  the  gravitating  body  itself.  Consider  a  gradiometer  which  can  be  mounted 
on  an  extendable  boom  so  that  measurements  can  be  taken  both  inside  that 
satellite  and  possibly  100-300  feet  away  in  the  direction  of  the  local  vertical. 
This  yields  the  gradient  of  the  gradient,  a  very  small  difference  in  output. 
Nevertheless,  it  is  theoretically  possible  to  determine  the  gravity  field 
Strength  from  these  two  measurements.  The  outputs  are  given  by 


Ag;  >  K  7T-  >/**<*" 


3.4- 22 

3.4- 23 


The  difference  between  these  two  outputs  is 

Sa?  -  K/l  (7(7  ~  yry) 


3.  4-24 
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This  may  be  approximated  closely  by 

&  T* 

$Af  «  3t C>  ~pr 


3.4-25 


where 


Gradiometer  scale  factor 

Nominal  distance  to  the  center  of  the  gravitational 
field 


dr  =  Boom  length 

Taking  the  sum  of  the  two  gradiometer  outputs  yields 

Za?  «  (>*rrbz) 

’I  Y*  I 


3.4-26 


Then  eliminating  rj  from  equations  (3.4-25  and  3.4-26)  and  solving  for yt( 
we  get: 


* 


27  (Z.A 
/€  K 


(Sr-)* 


3.  4-27 


All  parameters  are  theoretically  available.  Note  that  it  is 
now  possible  to  obtain  r,  explicitly  from  either  of  equations  (3.4-25  or 
3.4-26)  which  means  that  the  satellite  orbit  parameters  can  be  obtained. 

The  trouble  with  this  method  is  the  fact  that  Sa  g  is  so  small  that  small  errors  in 
A  g  give  large  errors  in  JA  g  and  any  errors  in  obtaining  that  latter  value 
are  further  magnified  by  the  sensitivity  factor  3  which  is  related  to  the 
exponent  of  g. 
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3.  4.  4  Torsional  Gradiometer 

An  instrument  which  relies  on  the  winding  or  unwind¬ 
ing  of  a  spring  as  compression  or  tension  is  applied  along  its  length  is 
presently  being  developed  by  the  Texas  Instrument  Corporation.  In  the 
vertical  orientation  tension  arises  as  the  difference  in  gravitational  attraction 
of  two  similar  proof  masses  mounted  at  either  end  of  the  spring.  The  con¬ 
sequent  angular  rotation  of  one  mass  (a  flat  plate)  relative  to  the  other  is  the 
measure  of  the  gradient. 

As  in  the  case  of  any  gradiometer,  the  general  advan¬ 
tages  of  this  instrument  is  its  insensitivity  to  externally  applied  acceleration, 
since  the  acceleration  acts  equally  on  each  proof  mass.  The  uncertain  future 
of  this  instrument  as  evaluated  by  Texas  Instrument  arises  from  its  early 
state  of  development.  Their  problems  include  size,  weight,  fragility  and 
stability  among  others. 

3.  4.  5  San  Marco  System 

The  San  Marco  experiment  concerns  the  observation 
of  the  motion  of  an  unsupported  object  within  the  confines  of  an  orbiting  body. 
The  motion  of  this  object  relative  to  the  body  is  due  to  the  fact  that  it  is  in  a 
different  orbit  from  that  of  the  body.  The  orbit  of  the  latter  can  be  defined 
in  terms  of  the  motion  of  its  center  of  gravity  relative  to  the  gravity  field. 
Since  the  test  object  is  not  located  at  the  center  of  gravity  its  orbit  differs 
from  that  of  the  satellite.  Thus,  the  object  will  move  within  the  satellite,  as 
a  free  object  in  space  until  it  strikes  an  obstruction  within  the  satellite.  By 
observing  the  motion  of  the  object  or  by  observing  the  magnitude  and  direction 
of  restraints  required  to  hold  it  fixed  relative  to  the  satellite,  the  gravity 
gradient  can  be  obtained. 

Consider  a  test  mass  attached  to  a  satellite  in  orbit  in 
a  central  force  field.  Let  the  test  object  be  alt  some  radial  distance  from 
the  satellite  c.  g.  Since  it  is  attached,  its  initial  velocity  components  are 
identical  with  those  of  the  satellite  (assume  no  rotation  of  the  satellite  about 
its  own  c.  g.  ).  It  should  be  noted  that  if  the  mass  were  sited  in  the  satellite 
on  the  trajectory  of  the  c.  g.  no  restraint  would  be  necessary  since  both  would 
then  be  on  the  same  orbital  path. 

When  the  test  object  is  released,  and  there  is  a  radial 
displacement,  the  test  object  is  not  in  the  same  orbit  as  the  satellite  and 
will  begin  to  move  With  respect  to  the  satellite.  It  is  this  motion  which  is 
the  measure  for  the  San  Marco  experiment  and  its  interpretation  is  shown 
in  the  following  analysis.  The  analysis  takes  into  account  the  effect  of 
releasing  the  test  body  with  a  non-zero  velocity  relative  to  the  satellite. 
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I  The  equations  of  motion  of  the  test  object  in  a  coordinate 

|  system  whose  origin  is  at  the  satellite  c.  g.  ,  rotating  at  the  angular  velocity 
w  of  the  satellite  are: 


- - 2*  -  xw2  ~  $ x/r 


3.  4-28a 


^+^x  +  xw-^wa  ~zfy/r 


3.  4-28b 


Here,  the  forces  producing  relative  motion  are  those  due  to 
the  displacement  from  the  c.  g.  multiplied  by  the  gravity  gradient  at  that 
point.  The  horizontal  gradient  is  one  of  compression  due  to  the  curvature 
of  the  equipotential  and  the  radial  gradient  is  one  of  tension  due  to  the 
distance  squared  effect  of  gravitational  force. 


Assuming  a  circular  orbit,  we  have  two  additional  relations: 


3.  4-29a 


g  =  -w2r 


3.  4-29b 


i  so  that  the  equations  of  motion  become 


x  -  2  wy  =  0 


3.  4- 30a 


y  +  2  wx  =  -  3y  w^ 


3.  4-  30b 


The  first  of  these  equations  can  be  integrated  using, 
as  initial  conditions,  x  =  xo  ,  Y  ~  Yo  >  x  =  u0,  y  =  .v0  at  time  zero. 


x  -  2  wy  =  u0  -  2  wyQ 


3.  4-31 


’  Sol  ving  for  x  and  substituting  in  (3.  4- 30b) 

##  o 

y  +  w^y  =  2  w  (2  wyQ  -  uq) 


3.  4-32 


This  undamped  second  order  system  cam  be  solved  by 
I  assuming  the  answer  is  in  the  form 


y  =  A  cos  wt  +  B  sin  wt  -  ?,  uq/w  +  4  yQ 


3.  4- 33a 


and  solving  for  the  coefficients  in  terms  of  the  initial  conditions. 
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Substituting  this  result  in  the  equation  (3.  4-31)  and  inte¬ 
grating  gives  the  horizontal  displacement  versus  time. 

x  -  x«  *  (+  Ut/w-6jj0)  sin  wt  +  2  CVt>/w)(l  -  oos  wt) 

,  #»  4  ,iu  \  3.4-33b 

-Jt 

It  should  be  noted  that  this  last  result  is  divergent  when, 
initially,  there  is  either  a  vertical  displacement  or  horizontal  velocity. 

For  the  test  where  only  a  vertical  displacement  exists  at 
time  zero  the  restilting  motion,  relative  to  the  satellite  eg  as  the  origin  of 
coordinates,  is  given  by 


6  yo  {  -  sin  wt  +  wt  ^ 
Yo  {  1  +  3  (1  -  cos  wt) 


3. 4- 34a 


3. 4 -34b 


which  is  a  spiraling  motion  in  the  x  direction  and  eventually  going  beyond  the 
bounds  of  the  physical  area  of  the  test. 


A  special  experiment  could  be  made  such  that 
u©  =  2  wyD 


3.4-35 


In  this  case,  the  divergent  term  in  x  would  be  zero  and  the  motion  of  the 
test  object  would  be  bounded.  The  resulting  elliptic  motion  is  given  by 

x  =  -2  y0  sin  wt  3. 4- 36a 

y  =  -yo  cos  wt  3.4-36b 

In  either  case,  gravity  gradient  is  determined  from  the  angular 
velocity  of  the  test  object  about  the  satellite  c.  g. 

Problems  associated  with  the  San  Marco  experiment  include: 

a)  the  problem  of  determining  the  center  of  gravity  of 
the  satellite 


b)  the  effects  of  the  change  in  satellite  eg  due  to  movement 

of  men  and/or  expendables  in  the  satellite. 

♦ 

c)  the  effects  of  solar  pressure,  micrometeorites,  etc. 
on  the  satellite  which  is  one  of  the  proof  masses  of 
the  experiment. 
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3.  4.  6  Gas  Diffusion  Technique  for  the  Determination  of 
Gravity  Gradient 

A  column  of  gas,  carried  on  board  the  orbiting 
satellite  will  take  on  a  density  gradient  related  to  the  gravity  gradient. 
Theoretically,  a  measure  of  that  density  variation,  probably  by  optical 
means,  will  yield  the  desired  data  on  the  gravity  gradient. 

The  force  exerted  by  the  gas  is 


F  =  pr  Vf>  3.4-37 

where  /o  (density)  and  p  (pressure)  are  related  by 
the  equation  of  state  for  an  ideal  gas. 

f-/>  £  T  3-4-38 

where  R  is  the  universal  gas  constant,  M,  the  average 
molecular  mass  and  T  is  the  absolute  temperature  of  the  gas.  Then 


%  rv  (iqf) 


3.4-39 


For  the  vertical  direction,  the  force  component  is 


Fi  * 


A  •  ol 


(  Ps 


3.  4-40 


3.  4-41 


This  is  to  be  balanced  against  the  force  of  gravity  in  the 

same  direction. 

e- ,  -  12  2  =  Z3°  r£  s  =  FT.  (Imp) 

2  °  (re+$)3  $  H  (  3 

This  can  be  integrated  (Peirce  #31)  to  obtain  the 
variation  of  pressure  as  a  function  of  altitude.  The  result  is 

>o 

In  exponential  form  we  have 

pfy  «  [- (m/ 3.  4-42b 

w 

Replacing  the  pressures  by  densities  from  the  equation 

of  state  one  obtains 

S  =/,o  *xf[-C*/*TXf?2/r»)l 


3.  4-42a 


3.  4-43 
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Thus,  by  a  careful  determination  of  the  variation  in 
density  of  a  suitable  gas  constrained  in  a  vertical  column  one  can  find  the 
variation  of  gravity  along  that  column. 

There  exists  a  variety  of  techniques  to  measure  the 
density  of  the  gas.  In  any  case,  the  detection  equipment  to  perform  such 
measurements  are  either  relatively  large  or  extremely  delicate.  Extreme 
care  must  be  taken  to  maintain  temperature  control  everywhere  in  the  gas 
column.  Temperature  variation  here  constitutes  the  major  source  of  "noise" 
which  may  mask  the  gradient. 

3.  4.  7  Summary  of  Gravity  Gradient  Instrumentation 

The  preceding  discussion  concerns  several  instruments 
and  systems  from  which  the  gravity  gradient  can  be  obtained  in  an  orbiting 
vehicle.  Instruments  include  the  Anna  Vibrating  String  Gradiometer,  a 
matched  pair  of  accelerometers,  the  gradiometer  on  a  boom  and  the  Texas 
Instruments'  Torsional  Gradiometer.  Systems  yielding  the  gravity  gradient 
include  the  San  Marco  system  and  the  Gas  diffusion  System. 

The  instrumentation  recommended  for  gradiometry 
measurements  is  the  vibrating  string  gradiometer  discussed  in  section  3.4.  1 
above.  This  recommendation  is  based  on  extensive  development,  manufactur¬ 
ing  and  testing  experience  with  the  closely  related  vibrating  string  acceler  - 
ometer.  Our  experience  leads  to  thg  conclusion  that  the  three  string  gradio¬ 
meter  can  be  produced  as  a  rugged,  compact,  highly  sensitive  device  for  the 
measurement  of  gravity  gradient  in  an  orbiting  satellite. 

An  error  analysis  of  the  gravity  gradient  system  shows 
that  the  vibrating  string  gradiometer  is  suitable  for  the  measurement  of 
gravitational  anomalies  in  the  order  of  ten  milligals.  A- growth  potential 
exists  which  will  substantially  improve  the  data  to  be  measured  from  an  orbit¬ 
ing  vehicle. 

3.  5  Systems  Analysis 

This  section  deals  with  the  manner  in  which  the  gradiometer  is 
introduced  into  the  lunar  gravity  environment  and  the  interrelationships 
among  the  instrument,  its  carrier  and  the  environments  which  affects  the 
determination  of  the  lunar  gravity  field. 

Specifically,  it  is  shown  that  relative  and  non-dimensional  data 
will  be  collected  and  that  orbit  parameters  may  be  obtained  from  the  instru¬ 
ment  in  this  environment.  An  example  of  a  gravity  measuring  system  is 
given  and  an  error  analysis  is  described  wherein  the  influences  of  the  system 
elements  on  the  system  accuracy  are  derived.  An  error  budget  gives  the 
results  in  quantitative  form  for  central  force  systems. 
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3.  5.  1  Relative  Versus  Absolute  Gravitational  Data 

An  important  aspect  of  gravitational  phenomena  is  the 
point  that  mass  is  relative,  and  its  magnitude  depends  on  the  value  chosen 
for  the  universal  gravitation  constant  G.  The  second  point  is  that  gravity 
gradient  technology  yields  data  which  have  the  dimension  of  time  and  we 
combine  them  with  an  independent  geometric  measurement  to  determine  the 
gravity  field  strength  and  related  functions. 


Consider  the  value  of  the  universal  gravitational  constant 
G  which  is  the  proportionality  constant  of  Newton's  Law  of  Mass  Attraction 


=  g  Ml  M2 

“72— 


3.  5-1 


It  may  be  worth  pointing  out  that  centuries  of  investi¬ 
gation  have  failed  to  yield  a  precise  value  of  G  (i.  e.  precision  to  within  a 
few  parts  per  million).  Instead,  the  current  value  of  G  is  "adopted". 
Assuming  distance,  r,  can  be  measured  to  the  precision  required,  the 
adoption  of  a  value  of  G  forces  the  magnitude  of  mass  to  fit  equation  (3.  5-1). 
In  other  words,  the  absolute  magnitude  of  mass  is  not  known  any  better  than 
\  the  absolute  value  of  G.  The  characteristics  which  can  be  found  to  a  high 
degree  of  precision  are  the  ratios  of  masses  and  the  functions  of  gravita¬ 
tional  field  strength  GM  namely  potential,  attractive  force  and  gradient, 
i  The  ability  to  determine  data  on  the  product  GM  results  from  the  fact  that 
tests  are  made  in  the  precisely  determined  dimensions  of  time  and  length. 

'  Thus,  two  masses  can  be  compared  by  "weighing”  them  in  the  same 

(Earth's)  gravitational  field  without  knowing  the  magnitude  of  that  field. 
Nothing  will  be  determined  in  absolute  terms.  Gravitational  force  is  deter¬ 
mined  by  noting  the  acceleration  to  which  a  faUing  test  body  is  subjected. 
This  is  independent  of  the  mass  of  the  test  body  as  Gaiileo  showed  long  ago. 

One  method  for  determining  the  ratio  of  the  masses  of 
the  Moon  and  the  Earth  has  been  the  determination  of  the  location  of  the 
center  of  gravity  of  the  Earth-Moon  system  which  is  in  orbit  around  the  Sun. 
It  is  not  the  center  of  the  Earth  itself  which  follows  an  elliptic  path  about  the 
Sun  but  rather  the  center  of  gravity  (barycenter)  of  the  Earth-Moon  system. 
Both  the  Earth's  center  of  gravity  and  the  Moon's  center  of  gravity  revolve 
about  the  barycenter  once  per  month.  From  astronomical  observations  of 
displacements  of  celestial  bodies  periodically,  relative  to  the  orbital  motion 
of  the  Earth,  the  barycenter  is  found  to  be  at  an  average  distance  of  2903 
:  miles  from  the  center  of  the  Earth.  The  eg  of  the  Moon  is  81.  3  times 

further  from  the  barycenter  Hence  the  ratio  of  the  masses  is  also  1:81.  3. 
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The  ratio  of  the  masses  can  also  be  found  from  the  con¬ 
stants  of  precession  and  nutation  of  the  Moon  in  the  Earth's  gravitational 
field.  H.  S.  Jones,  in  1941,  found  the  value  of  1:81.27  i0.02.  Then,  the 
accuracy  to  which  the  Moon's  mass  is  known  is  equal  to  that  to  which  the 
Earth's  mass  is  known.  Similarly,  one  can  determine  the  magnitude  of 
the  "average''  lunar  gravitational  acceleration. 

Obviously,  the  best  way  to  measure  the  Moon's  gravi¬ 
tational  field  is  by  use  of  gravimeters  located  at  various  points  of  the  Moon's 
surface,  just  as  it  isdone  on  Earth.  However,  this  is  not  feasible  at  the 
present  time.  The  next  best  method  would  involve  tests  from  a  lunar  orbit¬ 
ing  satellite. 

The  present  state-of-the-art  of  satellite  orbiting  is 
just  about  at  the  stage  where  a  vehicle  could  be  put  into  a  reasonably  precise 
orbit  around  the  Moon.  Such  an  experimental  station  would  greatly  enhance 
man's  knowledge  of  the  Moon's  characteristics  and  it  is  a  mission  that  should 
be  performed  before  landing  a  man  on  the  Moon.  Many  experiments  would 
be  conducted  from  an  orbiting  satellite  including  the  determination  of  atmos¬ 
phere  content  (to  the  extent  to  which  it  exists)  radiation  levels,  magnetic  and 
gravitation  fields.  Complete  three  dimensional  mapping  would  be  undertaken 
as  well  as  absorption  and  reflection  experiments  to  determine  the  nature  of 
the  moon's  surface. 

It  is  possible  to  sense  the  strength  of  the  gravitational 
field  from  within  an  orbiting  satellite  with  accelerometers  sensing  the 
gradient  of  the  field.  These  measurements  correlated  with  other  appropriate 
measures  will  yield  the  gravitational  field  strength  of  the  Moon  without  the 
necessity  for  landing  on  it. 

A  problem  with  regard  to  the  Moon's  gravitational  field 
is  its  non- spherical  shape.  This  causes  the  gravitational  field  to  be  a 
function  of  angle  as  well  as  radius.  The  fact  that  the  Moon  shows  approx¬ 
imately  the  same  face  to  the  Earth  throughout  the  month  indicates  that  it  is 
subject  to  gravity  gradient  attitude  stabilization.  That  is,  where  a  satellite 
has  unequal  moments  of  inertia,  the  gravity  gradient  effect  will  cause  the 
axis  of  least  moment  of  inertia  to  be  aligned  toward  the  center  of  the  gravi¬ 
tating  body,  i.  e.  along  the  local  vertical.  It  is  possible  to  improve  our 
knowledge  of  the  shape  of  the  Moon  by  gravity  field  measurements.  Jeffreys, 
1937,  has  shown  that  the  longest  axis  of  the  Moon  is  3478.  43  Km,  the  polar 
diameter  is  3-76.  25  Km  and  the  equatorial  diameter  normal  to  the  longest 
axis  is  3475.  55  Km.  Equipotential  surfaces  will  follow  this  contour.  Note 
that  the  Earth's  corresponding  dimensions  are:  long  axis,  in  equatorial  plane 
12,756.94  Km;  shorter  axis  in  equatorial  plane  12,756.  62  Km;  and  the  polar 
diameter  is  12,713.  82  Km.  Thus  the  departure  from  sphericity  of  the  Earth's 
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gravitational  field  is  normal  to  its  equatorial  plane.  It  does  not  have  a 
dumbell  analog  like  the  Moon.  The  Earth's  analog  is  a  toroid. 

3.  5.  2  Effect  of  Orbit  Eccentricity  on  Gradiometer  Output 


The  eccentric,  elliptic,  satellite  orbit  is  characterized 
by  a  variable  distance  between  the  satellite  and  the  center  of  the  gravity  field. 
The  locus  of  the  satellite  is 

r  =  A/(l  +  e  cos  0)  3.  5-2 

where  A  =  semi-latus  rectum 

0  =  true  anomaly 

e  =  eccentricity 

The  maximum  and  minimum  radial  distances  are  apogee 
and  perigee  (designation  with  respect  to  earth) 

ra  =  A/(  1  -  e)  (0  =  180°)  3.  5-3a 

rp  =  A/ ( 1  +  e)  (0  =  0o)  3.  5-3b 

The  radial  gradiometer  outputs  vary,  then,  during 
the  orbit  between  the  values 


Aga  =  K  (1  -  e)3 
A3 


3.  5-4a 


‘gp  -  K  ^3  (1  +  e)3 


3, 5-4b 


where  K  is  the  scale  factor. 


Therefore,  one  can  determine  the  eccentricity  of  the  orbit  from  the  cube 
root  of  the  ratio  of  maximum  to  minimum  outputs. 


(a^/av),/3 -0-'/(! 


■be)  ^  (j-e)1  ~  f-2 


3.  5-5 


for  small  eccentricities 


3.  5-6 
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This  is  a  non-dimensional  characteristic  and  can  be 
determined  from  gradiometer  measurements  alone.  Where  compensation  '■ 
for  eccentricity  is  required,  this  parameter  can  be  readily  determined. 


3.  5.  3  Example  of  System  for  Determining  Gravitational 

-Accelerate 


The  output  of  the  gradiometer  is  a  function  of  gravi¬ 
tational  field  strength  and  distance  from  the  center  of  the  gravity  field  in 
.  a  form  whose  dimension  is  the  reciprocal  of  time  squared.  In  order  to 
*  obtain  gravitational  acceleration,  the  distance  factor  is  introduced  explicitly. 
The  system  also  includes  a  gimbal  support  system  which  aligns  the  sensitive 
axis  of  the  gradiometer  to  the  local  vertical. 


Vertical  Error 


Orbit  }  r 


9  /  r~  \t 

Attitude 

G  radio  meter- 

0  f  1 

A, _ < 

Controller 

— nz 

U  c 

For  an  approximately  circular  orbit  the  gradiometer 

output  is  of  the  form 

^g  =  K0  +  Ki  £  cos2  C<  3.5-7 

Where  Kq  =  bias  of  instrument 

=  scale  factor  of  instrument 

^  =  attitude  of  sensitive  axis  relative  to 
local  vertical 


Assuming  that  the  orbiting  satellite  is  controlled  in 
attitude  to  the  local  vertical,  so  that  0(.  *  0  and  assuming  adequate  calibra¬ 
tion  methods,  the  gravity  field  can  be  measured  continually  in  orbit  s  a 
function  of  position  obtained  from  the  navigation  equipment. 


Thus,  a  map  of  the  gravitational  field  will  be  obtained, 
including  measures  of  anomalies  in  that  field  and  the  perturbations  due  to 
other  bodies  in  space. 
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As  a  first  order  model,  one  may  consider  that  the 
satellite  continually  reproduces  its  path  in  inertial  space.  It  is  then 
reasonable  and  easy  to  correlate  data  from  successive  circuits.  This 
correlation  will  improve  the  theoretical  accuracy  of  the  system  in  pro¬ 
portion  to  the  square  root  of  the  number  of  correlated  independent  data 
points. 

A  more  elaborate  model  such  as  suggested  in  section 
(3.  2.  2)  might  involve  a  three  dimensional  description  of  the  gravity  field 
in  terms  appropriate  to  the  least  squares  estimation  of  coefficients  in  the 
model.  This  would  require  more  tracking  data  and  would  eliminate  time  as 
an  explicit  parameter.  This  technique  is  similar  to  current  techniques  for 
evaluating  the  earth's  gravity  field  by  observation  of  satellites. 

In  addition  to  the  vertically  oriented  gradiometer,  the 
system  may  also  include  gradiometer s  aligned  normal  to  the  local  vertical 
in  the  orbit  plane.  The  output  then  is 

|  Lc*s*x-$™2*  sin1  0)  3.5.8 

This  output  indicates  the  orientation  of  the  gravity  equipotential ,  relative 
to  the  local  horizontal  and  the  angle  <p  from  the  orbit  plane.  For  a  homo¬ 
geneous  spherical  body,  the  equipotential s  would  be  concentric  spheres  and 
everywhere  coincident  with  the  geocentric  horizontal.  For  a  homogeneous 
oblate  or  prolate  spheroid  the  equipotential s  would  be  symmetric  to  the 
gravitating  body  near  its  surface  and  gradually  become  spherical  as  distance 
from  the  body  increased.  Where  inhomogeneities  in  the  figure  or  the  mass 
di  stribution  occur,  the  equipotential  will  be  most  distorted  near  those 
departures  from  symmetry,  gradually  disappearing  as  distance  increases. 
Thus,  for  a  gradiometer  maintained  aligned  to  the  geocentric  horizontal, 
in  a  satellite  in  orbit  close  to  the  surface  of  the  body  being  investigated, 
anomalies  will  be  detected  and  measured  for  relative  magnitude  and  extent, 
in  an  otherwise  noiseless  configuration.  By  holding  the  gradiometer  in  the 
nominally  null  position,  its  sensitivity  to  the  change  in  curvature  of  the 
equipotential s  is  very  high.  Note  that  the  third  gradiometer  aligned  normal 
to  the  two  mentioned  above,  will  give  the  third  dimension  to  the  magnitude 
and  orientation  measurement  of  the  equipotential  surfaces  through  which  the 
orbiting  satellite  passes.  This  system  will  yield  sufficient  data  to  test  the 
gravity  models  of  the  geophysicist. 
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3.5.4  Error  Analysis 


The  measurement  of  the  gravity  gradient  by  an  instru¬ 
ment  whose  attitude  is  maintained  along  the  local  vertical,  while  in  a  near 
ci  rcular  orbit,  is  combined  with  position  information  to  yield  gravity 
acceleration.  The  errors  in  such  a  system  are  determined  from  the  deriv¬ 
ation  of  the  system  equation  3.  5-7. 

S  -  r(A*-K,) 

o  yy  - -^Z  (Y 


K, 


COS' 


df  c  ft.  [(A}-K'e)‘Jr  +  r(‘jAf  d*,)] 
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3.  5-9 


3.  5-10 


3.  5-11 


CO s  oc 


r  '  *, 

The  functional  form  in  which  these  errors  appear  show 
the  weights  of  the  component  errors  on  the  desired  measure  for  the  system 
in  which  these  elements  are' assembled.  The  weights  associated  with  the 
errors  are  shown  parenthetically  in  the  following  five  relations. 


m(%)  3r 
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3. 5- 12b 
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small  angles. 
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Each  error  is  assumed  to  be  random,  zero  mean  and 
st ati stic ally  independent.  Therefore,  the  statistical  combination  of  the 
five  errors  given  in  the  above  equation  (3.  5-12)  is  a  vector  sum  rather 
than  an  algebraic  one. 


This  is  the  r.m.  s.  system  error  based  on  r.m.  s.  values 
of  the  individual  errors.  An  error  budget  can  be  established,  for  example, 
by  letting  each  error  source  contribute  equal  amounts  to  the  system  error. 
Equations  (3.  5-12)  show  that  unit  errors  in  the  individual  parameters  weigh 
differently  against  the  system.  Thus  the  budgeted  errors  should  be  in¬ 
versely  proportional  to  these  weights  in  order  that  the  parameter  errors  have 
the  same  effect  on  the  system. 

There  are  five  sources  of  error  noted  in  this  analysis. 

If  the  system  error  is  dg  and  the  contribution  of  each  source  is  dg,  the 
orthogonal  relationship  among  the  independent  error  sources  means  that 


dg  =  ypT  dg 


3.  5-14 


As  a  basis  for  numerical  evaluation  let  us  assume  that 
dg  should  be  determined  to  within  0.  02percent  of  its  true  value  (within  32 
milligals  on  the  Moon).  Then,  the  error  budgeted  to  each  error  source  is 
5  g  =  9  x  10-5g. 


Consider  that  the  measuring  system  is  in  a  nominally 
circular  orbit  about  100  miles  above  the  moon's  surface.  Then,  the  ex¬ 
pected  parameter  values  are 


g 

or 

Ki 

K0 

Ag 


=  1180  miles  (100  miles  altitude) 

=  136  gals  (at  altitude)  =  8.  5  x  10-4  mi/sec^ 

=  0  (attitude  controlled) 

=  1  (normalized) 

A  g  (by  design  and  calibration  techniques) 
=  Kj  g/r  =  7.2  x  10_7/sec2 
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Then,  from  the  nominal  values  and  budgeted  levels  the  tolerable  maximum 
errors  for  the  five  error  sources  are: 


*  r 

3  K0 

*Ki 


1180  x  9  10~5  _  mile 

7.2x10"  x9xlO“^  =  6.5xlO_ii/sec2 

6.5  x  10-11/ gec^ 


9  x  10-5 

-j/& 

(4.  5  x  10-5/ 


6.  7  x  10-3  radians  =  0.  38  degrees 


To  summarize  the  results  of  the  above  error  analysis, 
the  magnitudes  of  the  tolerable  errors  indicate  that: 


1)  The  gradiometer  should  be  so  designed  that  the  bias  is 
zero  theoretically. 


2)  Alignment  of  gradiometer  to  within  1/3°  of  the  local 
vertical  is  required. 

3)  The  accuracy  requirement  on  the  gradiometer  is  1/10,000 
in  this  example. 

It  has  been  assumed  that  the  gradiometer  should  be 
designed  for  a  linear  output  so  as  not  to  introduce  non-linear  error  sources. 
It  should  be  noted  that  by  taking  data  over  many  orbits  and  correlating  the 
results,  the  accuracy  of  the  desired  measurement  will  improve  as  the  square 
root  of  the  number  of  repetitions. 


3.  6  Design  of  the  VSG 

In  the  following  section  we  summarize  the  major  choices  and 
relationships  employed  in  the  design  of  a  Vibrating  String  Gradiometer,  as 
described  in  section  3.  4.  1.  Note  that  we  have  drawn  heavily  on  past  experience 
with  accelerometer  designs  for  BSD,  ASD,  and  NASA  to  assure  the  reason¬ 
ableness  of  the  selections.  It  is  concluded  that  a  VSG  can  be  made  to  meet 
the  requirements  of  lunar  anomaly  mapping  as  determined  earlier  in  this 
report. 
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3,  6.  1  String  Characteristics 

All  strings  are  made  of  the  same  material,  namely 
beryllium -copper.  The  density  of  this  metal  is  0.  3  lb/in3. 

The  strings  will  have  a  minimum  cross  sectional  area 
consistent  with  fabrication  and  handling  in  order  to  attain  the  highest  opera¬ 
ting  frequency.  Experience  indicates  A  =  0.  005  x  0.  001  square  inches  is 
reasonable. 


The  initial  design  for  the  end  strings  will  be  based  on 
Arma  experience  with  the  shortest  string  manufactured,  again  for  high 
nominal  operating  frequency  and  high  scale  factor  for  the  gradiometer. 

Lj  =  L>3  =  0.  375  inches. 

Further,  minimizing  the  end  strings,  minimizes  the 
overall  instrument  length. 


The  center  string  shall  be  as  long  as  possible,  to  obtain 
the  greatest  gradient  information.  Considering  handling  and  assembly 
problems  L<2  is  12  inches. 

The  nominal  operating  frequencies  of  the  three  strings 
are  approximately  3200  cps  for  the  end  strings  and  100  cps  for  the  center 
string. 


3.  6.  2  Proof  Mass  Characteristics 

The  proof  mass  is  made  as  large  as  possible  without 
overstre9sing  the  strings  or  the  cross  supports.  This  includes  the 
restriction  that  the  maximum  change  in  tension  M  (a  +  e)  is  small  compared 
to  T0  in  order  to  operate  in  the  linear  range  of  the  tension  versus  length 
relation. 


To  achieve  this  T0  in  the  selected  string  is  chosen  as 
10  grams  force.  This  is  the  least  T0  used  in  a  VSA  to  date.  Secondly  the 
proof  is  chosen  as  1.  2  kilograms  of  mass.  As  an  extreme  of  environment, 
leading  to  conservative  design,  the  tension  change  is  evaluated  for  the 
instrument  situated  1000  feet  from  the  zero-g  ambient  point,  the  center  of 
gravity  of  the  satellite,  at  1200  x  10~4  or  0.  1  gram  force  change.  Thus  the 
choice  of  the  mass  is  such  that  even  under  these  extremes  the  stress  is  well 
below  the  elastic  limits  of  Be  Cu  and  the  change  is  small  compared  with  the 
nominal  tension. 
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For  testing  purposes  in  the  Ig  environment,  however, 
a  separable  section  of  the  proof  mass  will  be  used.  This  test  mass  is 
about  6  gms.  and  will  be  unlocked  from  the  caged  proof  mass  for  the  test. 

The  matching  requirements  of  the  two  proof  masses  in 
the  gradiometer  is  important  in  keeping  the  bias  output  of  the  instrument 
well  below  the  level  of  the  nominal  output.  Theoretically  the  bias  is  zero. 
However,  in  the  presence  of  a  proof  mass  differential  dM  the  output  of 
the  unit  will  be  proportional  to 

Me  +  dM  (a  +  e) 

where  a  is  the  ambient  acceleration.  For  an  altitude  of  100  miles  above 
the  Moon,  for  example,  e  =1.4  x  10“ &  ft/sec^  for  an  instrument  one  foot 
long  and  a  =  1.4  x  10“ ^  ft/sec2  for  each  foot  radially  from  the  eg  of  the 
satellite  orbiting  the  Moon.  Assuming  the  gradiometer  were  100  feet 
radially  away  from  the  satellite  eg  it  is  necessary  that 

dM  «Me 
101  e 

Thus,  conservatively,  for  a  perturbation  of  less  than  0.  01  percent  dM  should 
be  less  than  M  x  10"  ^  or  0.  0012  grams.  The  masses  can  be  calibrated  by 
"weighing"  them  on  a  test  vibrating  string  since  VSA  experience  has  shown 
that  an  acceleration  (force)  resolution  of  this  order  of  magnitude  is  feasible. 

3.  6.  3  Cross  Coupling  Restraints 

With  regard  to  cross  coupling  of  energy  between  strings, 
the  structural  support  members  will  be  designed  with  internal  damping  and 
variable  cross  section  so  as  to  attenuate  any  energy  which  could  be  trans¬ 
mitted  from  one  end  string  to  the  other.  Such  cross  coupling  tends  to  reduce 
the  difference  frequency  of  the  two  strings  with  the  consequent  output  error. 

Cross  support  tapes  constrain  the  proof  masses  in  the 
two  directions  perpendicular  to  the  sensitive  axis  of  the  instrument.  Design 
analysis  shows  that  the  stiffness  of  the  proposed  cross  supports,  in  the  axial 
direction  will  be  only  one-three  hundredth  of  the  stiffness  of  the  center 
vibrating  string,  the  latter  being  Li/L2  as  stiff  as  the  end  strings.  Thus 
the  cross  supports  do  not  interfere  to  any  extent,  with  the  sensitivity  of  the 
instrument.  In  the  directions  normal  to  the  sensitive  axis  these  same  cross 
supports  are  25  x  103  times  as  stiff  and  will  prevent  change  in  the  orientation 
of  the  sensitive  axis  due  to  ambient  cross-acceleration  in  the  orbiting  vehicle. 
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3.6.4  Magnet  Assemblies 

The  VSG  requires  a  permanent  magnet  assembly  with 
each  string.  The  end- string  chosen  is  one  whose  length  is  equal  to  that  of 
an  Arma  VSA.  The  same  magnet  used  in  the  VSA  is  appropriate  for  the 
VSG.  Its  length  is  approximately  80  percent  of  the  string  length. 

For  the  center  string,  the  magnet  length  will  be  only  a 
small  portion  of  the  string  length.  This  is  feasible  since  very  little  energy 
is  required  to  drive  the  high  Q  strings.  The  magnet  assembly  for  the  center 
string  may  be  made  up  of  several  sections  placed  along  the  string  or  one 
magnet  at  the  center.  This  choice  will  be  optimized  during  the  final  design. 

3.  6.  5  Beaching 

The  instrument  is  protected  against  the  high  g  environ¬ 
ments  prior  to  being  set  in  orbit  by  the  use  of  automatic  beaching  (caging) 
of  the  full  proof  mass  at  a  point  below  the  elastic  limit  of  the  sensitive 
el  ements. 

3.  6.  6  Accuracy  of  the  VSG 

Generally  speaking  the  measured  output  of  the  VSG  is  a 
sum  of  frequencies  squared.  It  refers  to  a  bias  term  K0  an  acceleration  e 
(gradient  times  length)  times  a  scale  factor  Kj  and  a  residual  R  due  to  higher 
order  terms.  Analytically  it  is  of  the  form 

f2  =  K0  +  Ki  e  +  R  3.  6-1 

By  virtue  of  design  and  calibration  techniques  K0  and  R 
are  negligible  and  the  error  associated  with  the  determination  of  the  true  value 
of  e  is  due  to  the  errors  in  measured  f  and  calibration  of  Ki 

-15-  -£(*«  -¥--<2-r 
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3.  6-3 


The  instrument  error  is  the  square  root  of  the  sum  of 
squares  of  these  independent  errors.  Both  £f/f  and  ^Kl/Kj  are  obtainable 
in  the  order  of  10-5  or  better.  Then  the  resultant  uncertainty  is  .JsTx  10_5  = 
2.  2  x  10-5  This  can  be  reduced  by  longer  counting  periods. 


54 


DS  64-R531  -35 


3.6.7  Temperature  Sensitivity 


as  follows: 


The  main  effects  due  to  an  increase  in  temperature  are 


a.  The  outside  aluminum  structure  of  the  instrument 
increases  its  length  due  to  its  thermal  coefficient  of  expansion.  This  elong¬ 
ation  causes  an  increase  in  the  tape  tensions.  Increased  tape  tensions  cause 
increased  tape  frequencies. 

b.  The  inside  beryllium  copper  structure  increases 
its  length  due  to  its  thermal  coefficient  of  expansion.  This  expansion 
causes  a  decrease  in  the  tape  tensions.  Decreased  tape  tensions  cause  de¬ 
creased  tape  frequencies. 


By  arranging  these  effects  to  be  self  compensating, 
it  is  possible  to  design  the  instrument  to  be  relatively  insensitive  to 
steady  state  temperature  changes. 

Through  the  device  of  temperature  control,  and  through 
the  use  of  a  drift  model  in  the  data  processing  computations  these  effects 
can  be  further  reduced  to  admit  extended  intervals  between  recalibration 
operations. 


3.  6.  8  Summary  of  Performance  of  the  VSG 


Threshold  Output  10-13  g /ft 

Resolution  1/100,000  based  on  1  minute  counting; 

Low  noise. 

Range  10-  1  3  to  10- 6  g/ft 


Caging  Level 
Accuracy  (de/e) 
Size 
Weight 


10-4  g/ft 

2,  2  x  10-5  as  per  parai  3.  £,t 
Length  16  1/2",  Diameter  4  1/2" 
20.  5  lbs 
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3.  6.  9  Summary  of  Physical  Design 

In  summary  design  studies  and  past  experience  with 
similar  devices  show  that  the  proposed  VSG  is  feasible.  It  will  be  necessary 
to  extend  some  of  the  fabrication  techniques,  particularly  those  for  matching 
strings  and  masses  because  of  the  low  levels  of  acceleration  and  high 
accuracy  requirements. 

Since  temperature  variation  is  a  problem  in  any  sensi¬ 
tive  instrument,  it  may  be  desirable  to  employ  a  thermostatically  controlled 
constant  temperature  environment. 

The  computational  and  data  processing  functions  indicated 
for  the  VSG  output  are  but  a  small  part  of  the  computer  requirements  for  the 
geophysical  interpretation  of  the  data  so  collected.  It  is  recommended  that 
the  individual  string  frequencies  be  telemetered  to  an  Earth-based  station 
and  the  data  processed  at  a  convenient  computer  facility,  along  with  the 
geophysical  computations. 

3.  7  Instrument  Testing  and  Calibration 

The  vibrating  string  gradiometer  has  been  designed  to  operate  in 
a  very  low  ambient  acceleration  compared  to  the  one  g  environment  on  the 
ground.  It  is  necessary,  however,  to  test  the  unit  prior  to  actual  use  in 
orbit  to  give  assurance  that  the  instrument  is  in  operating  condition.  Once 
in  orbit  the  equipment  will  operate  as  planned.  Periodically  during  the 
course  of  the  satellite's  lifetime,  the  VSG  will  be  recalibrated  automatically 
for  bias  and  scale  factor.  This  will  enhance  the  accuracy  of  the  accumulated 
data. 


3.7.1  Prelaunch  Testing 

It  is  considered  appropriate  to  test  the  gradiometer 
prior  to  launch  to  the  fullest  extent  possible  even  though  the  ground  environ 
ment  is  recognized  to  be  far  beyond  the  extreme  limits  of  the  instrument. 
This  is  done  primarily  to  check  the  functional  operation  of  the  instrument. 

(Note:  It  is  possible  to  manufacture  the  equipment  in  a 
caged  configuration  and  release  the  proof  masses  only 
after  the  satellite  is  in  orbit.  All  testing  would  then 
be  done  on  the  satellite.) 
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The  gradiometer  can  be  operated  on  the  ground  in  a 
manner  analogous  to  its  operation  in  orbit  by  employing  test  masses. 

These  are  attached  to  the  strings  but  detached  from  the  caged  proof  masses 
during  test.  Design  studies  have  indicated  that  test  masses  of  about  one- 
two  hundredth  of  the  proof  mass  is  practical.  Then  tests  for  bias,  scale 
factor  linearity  and  cross  acceleration  effects  may  be  carried  out.  The 
gradiometer,  in  t«st  condition,  (proof  masses  caged,  test  masses  released) 
is  mounted  on  a  precision  rotary  head. 

When  the  VSG  is  mounted  normal  to  the  local  vertical, 
the  frequencies  of  the  two  end  strings  should  be  identical.  The  difference 
in  frequency  of  the  two  end  strings  is  compared  to  that  difference  when  the 
unit  is  rotated  180°  about  the  local  vertical,  i.  e.  with  the  sensitive  axis 
still  horizontal.  One  half  the  sum  of  these  differences  is  the  bias,  since 
the  rotation  takes  account  of  the  possible  initial  misalignment. 

Scale  factor  and  linearity  in  the  acceleration  measuring 
mode  are  determined  by  testing  the  output  of  the  instrument  at  several  points 
in  the  one  g  gravity  field.  Here,  the  gradiometer  is  rotated  about  the 
horizontal  normal  to  its  sensitive  axis.  The  sine  of  the  angle  from  horizon¬ 
tal  is  the  portion  of  the  gravity  field  being  measured. 

For  calibrating  the  gradiometer  function,  the  instru¬ 
ment  is  placed  on  a  horizontally  rotating  platform  (centrifuge)  and  the  output 
determined  as  a  function  of  angular  velocity  and  position. 

The  effect  of  cross  acceleration  on  the  gradiometer 
output  is  determined  by  noting  the  variation  in  output  as  the  gradiometer 
is  rotated  about  its  sensitive  axis  while  this  axis  is  held  horizontal.  The 
instrument's  cross  support  system  is  the  primary  mechanism  by  which 
cross  acceleration  effects  are  minimized. 

3.7.2  In-Orbit  Calibration 

The  prelaunch  tests  assure  that  the  instrument  and  its 
associated  elements  are  functionally  operative.  Since  the  actual  proof 
masses  are  not  used  in  the  one  g  environment,  calibration  is  required  in 
orbit.  There  are  several  possible  techniques  which  may  be  employed  to 
determine  the  bias  and  scale  factors  of  the  instrument. 

In  orbit,  the  test  masses  (which  are  secured  to  the 
sensitive  strings)  are  locked  to  the  proof  masses  and  the  combination  will 
be  uncaged.  Assuming  that  the  gradiometer  is  mounted  in  its  own  gimbal- 
ing  system  so  that  its  orientation  can  be  modified  automatically,  the 
gradiometer  will  be  rotated  about  its  own  eg  at  various  angular  rates.  The 
change  in  output  as  a  function  of  angular  rate  gives  the  scale. 
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The  bias  output  would  be  determined  in  a  manner 
analogous  to  that  outlined  for  the  prelaunch  test.  It  is  important  to  note 
that  the  bias  may  also  be  isolated  in  geophysical  data  processing. 

A  second  possible  method  for  calibrating  the  gradio- 
meter  in  orbit  employs  the  effect  of  a  known  mass,  on  board  the  satellite, 
on  the  gradiometer  output.  Consider  a  "calibration  mass"  suspended 
along  the  sensitive  axis  of  the  gradiometer  on  the  side  opposite  that  of  the 
earth  (i.  e.  gravitating  body).  The  presence  of  this  mass  in  close  proximity 
to  the  instrument  should  be  sufficient  to  offset  the  effect  of  the  earth's  gravity. 
The  distance  between  the  test  mass  and  instrument  when  the  output  of  the 
latter  is  zero  represents  the  scale  factor  (after  the  bias  has  been  taken  into 
account).  Further,  by  noting  the  outputs  as  a  function  of  distance  the  input- 
output  function  is  determined. 

There  is  a  monitoring  technique  which  can  be  used 
without  any  special  constraints  on  the  gradiometer  and  which  is  effective 
all  during  the  normal  operating  time  of  the  instrument.  This  relates  the 
maximum  and  minimum  outputs  of  the  unit  to  the  period  of  the  satellite  in 
an  eccentric  orbit. 
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3.  7-1 


3.  7-2 


The  time  interval  between  successive  perigee  passages,  which  can  be 
measured  inboard,  is  related  to  the  orbit  by 


T- 2ir  p/V 

where  a  =  semi-major  axis  =  1/2  (Tcl  +  70 
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The  relative  distances 
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are  thus 


related  by 

3.7-5 


This  latter  data  can  be  extracted  from  the  system  and  the  equality  checked 
each  orbit  period. 
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